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Abstract 

The spectral energy distributions (SEDs) of galaxies are shaped by nearly every physical 
property of the system, including the star formation history, metal content, abundance 
pattern, dust mass, grain size distribution, star-dust geometry, and interstellar radiation 
field. The principal goal of stellar population synthesis (SPS) is to extract these variables 
from observed SEDs. In this review I provide an overview of the SPS technique and discuss 
what can be reliably measured from galaxy SEDs. Topics include stellar masses, star 
formation rates and histories, metallicities and abundance patterns, dust properties, and 
the stellar initial mass function. 



1 INTRODUCTION 

Many of the fundamental properties of unresolved stellar populations are encoded in their 
SEDs. These properties include the star formation history (SFH), stellar metallicity and 
abundance pattern, stellar initial mass function (IMF), total mass in stars, and the physical 
state and quantity of dust and gas. Some of these properties are easier to measure than 
others, and each provides important clues regarding the formation and evolution of galaxies. 
It is precisely these quantities, measured from the SEDs of galaxies, that have provided the 
foundation for our modern understanding of galaxy formation and evolution. 

Over the past several decades considerable effort has been devoted to extracting infor- 
mation from the SEDs of galaxies, exploiting information from the FUV to the FIR. Early 
attempts at understanding the visible and NIR spectral windows approached the problem 
by co mbining mixtures of sta rs in ad hoc ways until a match was achieved with observations 
(e.g JSpinrad fc Tavloll97lh . More sophisticated versions of this techniq ue were devel oped 
that incorporated physical constraints and automated fitting techniques (|Faberl ll972). At 
about the same time, synthesis models were being developed that relied on stellar evolution 
theory to constrain the range of possible stella r types at a given age and metallicity (e.i 



tncory to constrain tnc range ot possible stellar types at a given age ana metallicity 
Tinslevl[l96sl . ISearle. Sargent fc Bagnuolclll973l . iTinslev fc Gunnlll97rj . lBruzuaJ[l983h 



The 



substantial progress made in stellar evolution theory in the 1980s and 1990s paved the way 
for t he latter approach to become the de facto standard in modeling, the SEDs of galaxies 
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modeling technique, which will be described in detail in the next sect ion, is sometimes re- 
ferred to as 'evolutionary population synthesis' (e.g.. lMarastonlll998T ). although the term 



'stellar population synthesis' (SPS) has garnered wider use. The latter term will be used 
throughout this review. 

The UV and IR spectral windows are rather more difficult to probe owing to the obscuring 
effects of the atmosphere. Nonetheless, numerous balloon and space-based observatories 
have opened up the ultraviolet and infrared to detailed investigations. In these spectral 
regions dust plays a major role; it absorbs and scatters much of the UV light emitted 
by stars and re-radiates that energy in the IR. In young stellar populations the UV is 
dominated by hot massive stars, while in old stellar populations the UV can be influenced 
by hot evolved stellar types such as post-AGB and extreme HB stars (see the review by 
0'Connelj[l999l for details). 

Th e development of models for the IR SEDs of galaxies (e.g., Draine fc Leejl984 . Zubko. Dwek fc Arendt 



2004) has proceeded in parallel with the development of models for UV, optical, and NIR 
SEDs, and it is only recently that models have been developed to simultaneously and self- 



consistently predict the FUV through FIR SEDs (e.g ., Devriendt. Guiderdoni fc Sadatll 999: 



Silva et al.lll998l ; Ida Cunha. Chariot fc Elbazlboosl ; iGroves et allboosl ; INoII et al.|[2009al ') 

Several broad questions will serve to focus this review. They are: 

• What have we learned about the physical properties of galaxies from their observed 
SEDs? 

• How reliable are the quantities thus derived? 

• What can be learned, in principle, from the SEDs of galaxies? 

The first question is relatively straightforward to address, while the second and third are 
necessarily more complicated. Attention will be given to cases where answers to these 
questions are presently unknown, but knowable (the 'unknown unknowns', in Rumsfeld's 
sense, are of course the most interesting, but most difficult to discuss). 

There are surprisingly few thorough reviews of stell ar population s ynthesis and its applica- 
tion to modeling SEDs. The foundational review bv lTinslev |l980h is highly recommended 
to anyone seeking a broad yet in t uitive und erstand i ng of stellar populations and galaxy 
formation. The reviews bv lFaberl (|l977r ) and lFrogell (|l988l ) are more narrowly focused on 
old stellar populations but they a re informative becaus e many of the issues raised therein 
are still relevant today. Recently, IWalcher et al] (|201lh has presented an excellent review 
of many aspects related to the modeling of galaxy SEDs. 

The topic of modeling galaxy SEDs is vast, and it would be impossible to provide a 
thorough review of the entire field. Hard decisions therefore had to be made. With regard 
to wavelength, energies higher than the FUV (roughly the Lyman limit), and lower than 
the submillimeter (roughly 1mm) will not be discussed. Nebular emission lines will also 
be neglected, except in a few cases. The measurement of and uncertainty induced by 
photometric redshifts will not be discussed and active galactic nu clei (AGN) will be ign ored. 
The technique of fitting models to data was reviewed recently bv lWalcher et al] (|201 In and 
so will not be discussed here. A whole review could (and should) be written on the evaluation 
and comparison of existing SPS models. This will not be undertaken herein, except in cases 
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where differe nt models, when applied to data, produce starkly different results. The reader 
is referred to Conrov fc Gunn 1 201ph for a recent comparison of several popular models. 
Finally, it is worth emphasizing what this review is not: it is not a summary of science 
results derived with SPS models. There is an extensive literature devoted to applying SPS 
models to data in order to derive insights into the formation and evolution of galaxies. 
Results of this nature, though fascina ting in their ow n ri ght, will not feature promine ntly 
in this review. Readers are referred to Renzini 1 20061 ) and Blanton fc Moustakas ( 20091 ) for 
recent reviews along these lines. 

This review is organized as follows. An overview of SPS model construction and ap- 
plication is provided in Section 2. We then turn to the topics of mass-to-light ratios and 
stellar masses (Section 3), star formation rates, histories, and stellar ages (Section 4), stellar 
metallicities and abundances patterns (Section 5), dust (Section 6), and the initial mass 
function (Section 7). Several concluding remarks are provided in Section 8. 



2 OVERVIEW OF STELLAR POPULATION SYNTHESIS 

The construction of models for simple and composite stellar populations is conceptually 
straightforward. There are however certain constraints that make the creation of such 
models rather difficult in practice (incomplete isochrone tables, incomplete empirical stellar 
libraries, poorly calibrated physics, etc.). In this section the ingredients necessary for 
constructing model SEDs will be discussed. Areas of particularly large uncertainties will be 
highlighted, but it is not the goal of this review to provide an exhaustive inter-comparison of 
various possible model choices. An overview of the entire process of constructing composite 
stellar populations is given in Figure [1] 



2.1 The Simple Stellar Population 

The starting point of any SPS model is the simple stellar population (SSP), which describes 
the evolution in time of the SED of a single, coeval stellar population at a single metallicity 
and abundance pattern. An SSP therefore requires three basic inputs: stellar evolution 
theory in the form of isochrones, stellar spectral libraries, and an IMF, each of which 
may in principle be a function of metallicity and/or elemental abundance pattern. These 
components are typically combined in the following way: 



fssp(t,Z) 



,(t) 



/star [T cB (M), log g(M)\t, Z] $(M) dM, 



(1) 



where M is the initial (zero-age main sequence) stellar mass, $(M) is the initial mass 
function, / s t ar is a stellar spectrum, and /ssp is the resulting time and metallicity-dependent 
SSP spectrum. The lower limit of integration, mi , is typically taken to by the hydrogen 
burning limit (either 0.08 or O.1A^0 depending on the SPS code), and the upper limit is 
dictated by stellar evolution. The isochrones determine the relation between T e g, logg, 
and M for a given t and Z. This approach to constructing SSPs is common but not 
universal. Alternatives include the fuel consumption theorem ( Renzini fc Buzzoni ,1986, 
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dBica & Alloinl 
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Ne arly all SPS mode ls provide SSPs as a 'black-box' (e.g., Leitherer et al .1 19991 . iBruzual fc Chariot 



20031 . lMarastonll2005l . IVazdekis etHlboioh . The user therefore has little working knowl 
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edge of how SSPs are built and where major issues lie. The following discussion is therefore 
largely pedagogical and geared toward users of SPS models, rather than model builders. 

2.1.1 Stellar Evolution & Isochrones An isochrone specifies the location in the 
Hertzsprung-Russell (HR) diagram of stars with a common age and metallicity. Isochrones 
are constructed from stellar evolution calculations for stars from the hydrogen burning limit 
(~ O.IMq) to the maximum stellar mass (« 100 M@). The construction of isochrones is 
straightforward for stellar evolution tracks that are infinitely well-sampled in mass and 
time. In practice, evolutionary tracks are discretely sampled and this can lead to issues 
in isochrone construction for fast evolutionary phases. Modern sets of isochrones have 
constructed specifically t o ensure that the models are relatively immune to these effects 



( Chariot fc Bruzuallll991 



A number of widely-used isochrone tables exist in the literature. The most popular mod 
els span a wide range in ages (masses), chemical compositions, and cover most relevant 
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Girardi et al. 


200C, 
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2007) 



2000l. iMarigo et all 120081 ) . an d the BaSTI models JPietrinferni et al 



20071 ). The Geneva models jSchaller et all 1 19921. iMevnet fc Maeder 



20001 ) are tailored to follow high-mass stars through advanced evolutionary phases, includ 
ing the Wolf-Rayet (WR) phase, but they do not model low mass stars. Other mod 
els have focused on the main sequence, red giant branch (RGB), and ho rizontal branch 
(HB) evolution of low-mas s stars (M < 3 Mo), incl uding the Y 2 mod els l Yi et al.l 12001 



Yi. Kim fc Demarquell2003l). the Dartmou t h models dDotter et alJboojf). and the Vi ctoria- 
Regina models ( Vandenberg fc Bell 1985 . VandenBerg. Bergbusch fc Dowler 20061 ) . Fi- 
nally, there are isochrones tailored to very low-m ass stars and brown dwarfs. The Lyon 
mode ls are the most widely used in this regime ( Chabrier fc Baraffei 1997 . Baraffe et al. 
19981 ). It is noteworthy that none of the isochrones listed here cover the post- AGB evolu - 



tionary phase. Remarkab ly, th e post- AGB iso chrones computed by Schoenberner ( 19831 ). 
Vassiliadis fc Wood! (1 19941 ). and iBloeckerl (|1995l ) are still widely used in modern SPS codes. 

Implementing isochrones in an SPS model is challenging because no single set spans the 
necessary range of ages, metallicities and evolutionary phases. It is common to use the 
Padova isochrones for the bulk of the age and metallicity range and to supplement with the 
Geneva models at young ages. Little attention is normally paid to the lowest mass portion 
of the isochrones, since low mass stars contribute on ly ~ 1% of the light of an old stellar 
population. An exception in this regard is the model of lConrov fc van Dokkuml (|2012al ). who 
paid special care to the modeling of low-mass stars for SPS. The splicing together of various 
isochrone sets can be difficult. For example, different codes make different assumptions 
regarding convection, rotation, etc., and so the age at which certain stars evolve off the 
main sequence varies between codes. 

Modules for Experiments in Stellar Astrophysics (MESA) is a new, highly modular and 
sophisticated stellar evolution code that includes the latest stellar interior ingredients in- 
cluding opa city tables, equatio ns of state, nuclear reaction networks, and surface boundary 
conditions ( Paxton et al. 201 lh . There is great hope that MESA will be employed to pro- 
duce high-quality isochrones over the full age and metallicity range and for all evolutionary 
phases. 

In addition to the practical difficulty of implementing current isochrone libraries, stellar 
evolution calculations contain a number of uncertainties relevant to SPS. For example, 
all of the stellar models mentioned above are based on one dimensional codes. As such, 
they require approximations to inherently three dimensional processes such as convection, 
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Figure 1: 

Overview of the stellar population synthesis technique. The upper panels highlight the ingredients 
necessary for constructing simple stellar populations (SSPs): an IMF, isochroncs for a range of 
ages and metallicitics, and stellar spectra spanning a range of T e ff , Lboli & n d metallicity. The 
middle panels highlight the ingredients necessary for constructing composite stellar populations 
(CSPs): star formation histories and chemical evolution, SSPs, and a model for dust attenuation 5 
and emission. The bottom row shows the final CSPs both before and after a dust model is applied. 



rotation, mass-loss, close binary interactions, and thermal pulses during AGB evolution 
These processes lead to major uncertainties in the isochrones that impac t t he resulting; SPS 
model predictions (e.g. Jcharlot. Worthev fc Bressanlll996l ; ICharlotlll996l ; Iy|2003| ; iLee et al 



20071 ; IConrov. Gunn fc Whitell2009r i 



The cores of stars more massive than M ~ 1.1 Mq are convective, as are the envelopes of 
evolved giants and low mass stars. Classically, the boundary between convective and radia- 
tive regions is specified by the Schwarzschild criterion. However, this criterion is effectively 
a requirement that the acceleration of a convective fluid element be zero. The fluid element 
will likely have a non-zero velocity as it crosses this boundary, and so some amount of 'over- 
shooting' is expected. This will result in a larger convective region than would be expected 
from the Schwarzschild criterion alone. With regards to the convective core, a wide body 
of observational evidence favors t he existence of a moderate amount of overshooting in the 
mass range 1.1 < M/M B < 7 Jstotherslll99ll. iNordstroem. Andersen fc Anderser]|l997l . 



VandenBerg fc Stetson! 12004 . iKeller fc Wood! 120061 ) . The amount of core overshooting re 



quired to fit the data results in a ~ 25% increase in the main sequence lifetime compared 
to models that do not include overshooting. Nearly all SPS models adopt isochrones with a 
modest amount of overshooting in t he convect i ve core, as s upported by observations. The 
exception to this is the SPS model of lMaraston 1 19981 . [2005J), which uses isochrones without 
core convective overshooting. The treatment of core overshooting has a noticeable effect on 
the color evolution of SSPs in t he age range of ~ 0.1 — 1 Gyr, by as much as ~ 0.1 mag 
(Yil booi IConrov fc Grmnlboiol ) and can therefore be an important source of systematic 
uncertainty when modeling SED s. 

As discussed in Cassisi the amount of overshooting in the convective envelopes of 

evolved giants is less constrained. This is unfortunate because the treatment of envelope 
convection affects, a mong other observables, the ratio of red to blue helium-burning giants 



19921 ). These stars can contribu te several tens of perce nt to the integrated 
light, depending on the star formation history ( Melbourne et al. 20121 ). and so uncertainty 



(Renzini et al 



in the amount of envelope overshooting should affect the integrated light predictions at a 
significant level. Overshooting in the convective envelope also affects the location of the 
luminosity bump in the RGB luminosity f unction. This fact was used bv lAlongi et al.l (|1991) 
and more recently bv lCassisi et al.l (|2011T ) to argue for a modest amount of overshooting in 
the convective envelope. 

The importance of stellar rotation h as been investigated by the Ge neva group, amongst 
others, over the past two decades (see Maeder fc Mevnet 2000l . 20121 . for a review). Per- 
haps most importantly, rotation increases the main sequence lifetimes (by ~ 25%) due to 
rotation-induced mixing bringing fresh fuel to the convective core. Rotation also lowers 
the effective surface gravity, lowers the opacity in the radiative envelope, increases the lu- 
minosity, and changes the ratio of red-to-blue supergiants. The mixing to the surface of 
H-bur ning products caused by r otation will also affect the number and type of Wolf-Rayet 



Vazquez et al. (2007) and Levesaue et al. ( 20121 ') investigated the impact of rotating 



stars. 

massive star models on integrated light properties by comparing to non-rotating models and 
found that the number of ionizing photons increased by as much as an order of magnitude 
and colors became bluer by as much as 0.1 — 1 mag, depending on wavelength and age. The 
signatures of WR stars as a function of time and metallicity in integrated spectra are also 
different between rotating and non-rotating models. 

Most massive stars are in binary systems, and there is some evidence that massive star 
binar ies preferentially have comparable masses ( Kobulnickv fc Fryer 2007 . Sana fc Evans! 
201 in . The interaction of close binary stars via mass-transfer and common envelope evolu- 
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tion will affect t he evolution of such stars and b ring; about further changes to their observ- 
able properties. Eldridee. Izzard fc Tout ( 2008h demonstrated that in many respects the 
effect s of binary star mod e ls wit hout rotation are similar to single star models with rota- 



Eldridge fe Stanwav argued that SPS models including binary star evolution 



tion. 

produced a better fit to UV spectra of star-forming galaxies. Regardless of the details, it 
is clear that where massive stars matter in galaxy SEDs, the effects of both rotation and 
binary evolution will play an important role. In fact, binary star evolution may also cre- 
ate blue straggler stars and extreme HB st ars, which would suggest that binary evol ution 
can affect older stellar populations as well (|Han et al.ll20ol . hooj . Izhang et al.j|2005h . No 
popular SPS model includes the effects of binary star evolution. 

The potential importance of thermally-pulsa t ing A GB (TP-AGB) stars in the context of 
SPS models was emphasized by Maraston et al. 1 20061 ). and has since become a controversial 
topic. This phase occurs for stars in the mass range « 1 < M/Mq < 8 (depending on 
metallicity) and is difficult to model for several reasons, including the fact that nuclear 
burning occurs in alternate hydrogen-rich and helium-rich shells. When the helium shell 



burns it does so explosively because of the thin shell instability (jSchwarzschild fc Harm 



1968), which gives rise to the thermal pulses. Mass-loss becomes catastrophic during this 
phase, thereby terminating the life of the star. Recently the Padova group has developed 
a new suite of isochrones with upd ated TP-AGB models ca l ibrated against obs ervations in 
the Large Magel l anic C loud fLMC lMarigo fc Girardll2007l . iMarigo et aljboosj ). However, 
Conrov fe Gmml (|201ol ) found that the updated Padova models failed to reproduce the 
colors of intermediate-age Magellanic Cloud star clusters, by as much as 0.5 mag in some 
cases. These authors provided recalibrated SPS models in which the wei ght given to TP 



AGB stars was reduced in order to match the LMC data. More recently. [Melbourne et al 



( 20121 ') analyzed resolved color-magnitude diagrams in nearby dwarf galaxies and concluded 
that these updated models produce twice as much luminosity in the TP-AGB phase as 
observed, roughl y inde pendent of the inferred mass fraction in young stars. In addition, 
Melbourne et al] (|2012h found that the latest Padova model predictions for the luminosity 
contributed by red core helium burning stars is a factor of two higher than observed in 
dwarf galaxies. This problem in the models may be related to the treatment of convection 
in the envelopes of evolved giants. 

Mass- loss is another critical parameter in stellar evolution models. At M < 8Mq, it 
determines when a star will end its life as a white dwarf, and how massive the white 
dwarf will be. At higher masses mass-loss can significantly alter the course of advanced 
evolutionary phases, especially at M > 40 Mq. In high mass stars the mass loss mechanism 
is thought to be via line-driven winds, while in lower mass st ars (< 8Mq) it is believed to 
be due to pulsation-induced dust-driven winds ( Willson 2000l ) . In any event, the mass-loss 
prescription is another free parameter in these models, and it is a critical one because it 
strongly affects the lifetimes of advanced (and luminous) evolutionary phases. For example, 
the lifetime of TP-AGB stars and th e number of thermal pu lses they undergo depend 
strongly on the mass-loss prescription (jVentura fc Marigollioiol '). 



In summary, the computation of isochrones for use in SPS depends on many uncertain 
aspects of stellar evolution including the treatment of convection, close binary evolution, 
rotation effects, and mass-loss, amongst many others. The importance of these uncertainties 
on derived SPS results will be highlighted throughout this review. 



2.1.2 Stellar Remnants Stars eventually die, usually leaving behind stellar rem- 
nants in the form of white dwarfs, neutron stars, or black holes (theory predicts that a cer- 
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tain class of very massive m etal-poor stars un dergoes pair-instability supernovae that leave 



behind no stellar remnant; Heeer et al. 20031 ). The relation between the initial, zero-age 



main sequence stellar mass and the final remnant mass is not well-constrained observation- 
ally, especially for massive stars. The initial-final mass relation for white dwarfs can be 
reasonabl y well constrained b y measuring white dwarf masses in open clusters with known 
ages (e.g., Kalirai et al.l 20Q gf). The initial-fina l mass relation is predicted to be a function 



of metallicity i|Marigoll200ll , Irleger et~al]|2003h . further complicating the situation. 



In SPS models stellar remnants are usually included in the total stellar mass budget, and 
their contribution can be significant. For example, if a certain amount of mass is formed 
into stars instantly, then after 13 Gyr only 60% of the mass remains in stars or stellar 
remnants; the other 40% has returned to the interstellar medium. Of the remaining stellar 
mass, 25% is locked in stellar remnants. Of the total remnant mass, 73% is comprised 
of white dwarfs, 7% is in neutron stars, a nd 20% i s in b lack holes. These numbers were 
computed ass uming solar metall i city a nd a Kroupal (2001) IM F, with the initial-final mass 

SPS model predicts similar 



relations from Renzini fc Ciotti ( 19931 ). The 



Marastor 



nl |2005h 



numbers (largely because the same initial-final mass relation was used). The point is that 
stellar remnants can make an important contribution to the total stellar mass of a system. 

2.1.3 Stellar Spectral Libraries Stellar spectral libraries are required to convert 
the outputs of stellar evolution calculations — surface gravities, g, and effective tempera- 
tures, T e ff — as a function of metallicity, Z, into observable SEDs. There is however no 
single spectral library, whether theoretical or empirical, that covers the entire range of pa- 
rameter space necessary for constructing SPS models. Stitching together various libraries, 
often of widely varying quality, is therefore necessary. Some modelers take the approach of 
constructing spectral libraries entirely from theoretical calculations, while others use purely 
empirical libraries. The benefits and drawbacks of these two approaches will be discussed 
in this section. 



Theoretical Libraries Theoretical libraries offer the great advantage of densely covering 
parameter space, including spectral resolution, and of producing spectra that are not subject 
to observational issues such as flux calibration and atmospheric absorption. The clear 
disadvantage is that the libraries are only as good as the input atomic and molecular 
parameters and the approximations made in the computation of the models, as discussed 
below. 

There are a number of decisions that must be made when computing synthetic spectral 
libraries including how to treat convection and the microturbulent velocity profile, and 
whether or not to model departures from LTE and plane-parallel geometries. Additional 
important limitations includ e the incompl ete and often inaccurate atomic and molecular 
line lists. As emphasized by iKuruczl (|201ll ). models still fail to reproduce all the observed 
features in ultra high-resolution spectra of the Sun due to incomplete atomic line lists. The 
situation is even more serious for cooler sta rs because the molecular line lis ts can carry 
fairly large uncertainties, in particular TiO (|Allard. Homeier fc Frevtadl201ll ). A related 



problem is that a significant fraction of the atomic and molecular lines are derived from 
theoretical calculations, rather than measured in the laboratory, and so they have uncertain 
strengths and central wavelengths. Unfortunately, these 'predicted lines' can be important 
for determining the overall SED shape of stars and hence of galaxies. Worse still, the 
partition function for many molecules is poorly known, so even the total abundance of 
particular molecules is uncertain. 
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The quality and state of the art of the theoretical models varies considerably as a function 
of T e ff . At the high temperature e nd, i.e., Wolf-Rayet and O-type m ain se quence stars, the 
state- of-the-art libraries are from ISmith. Norris fc Crowther (2002ft and lLanz fc Hubenvl 
( 20031 ). while for hot t er com pact objects (e.g., post-AGB stars) the most up-to-date mod- 
els are from iRaucbl hooj). By far the widest range in parameter space is covered by 

This librar y 

(mi), 



the BaSeL atlas ( Leieune. Cuisinier fc Buser 1997 . 19981 ; Westera et al 



2002ft. 



Fluks et al 



is comprised of a variety of th eoretical models from Kuru cz (1995), 
Bessell et all (| 19891 . Il99lh . and lAllard fc Hauschildtl (|l995l V The broadband SEDs of the 
stars in this library have been modified to agree with observed UBVRIJHKL color- 
temperature relations for individual stars. The BaSeL atlas is almost universally used 
in modern SPS codes, despite the fact that the input theoretical spectra are now almost 20 
years old. At the very coolest temperatures the MARCS and NEXTGEN/PHOENIX grids 
are the state-of-the-art because of their comprehensive molecular line lists an d treatment of 



spheri cal effects on the atmospheres and spectra of cool giants. Finally, the lAringer et al 



(2009) models fill out the low temperature end by providing carbon star spectra over a 
range of parameter space. 

The re are several modern spec tral libraries computed bv lMunari et al. I l|2005ft .l iMartins et al 
(2005J), and lCoelho et all ( 2005T ) that offer fairly wide coverage in parameter space and are 
at high spectral resolution. The IMartins et al.l (12005) library has even been in corporated 



into a fully theoretical high resolution SPS model (|Gonzalez Delgado et al . 2005). However 
these libraries appear to be geared toward high resolution spectroscopic analyses because 
they do not include the predicted lines. They therefore cannot be used to model broadband 
SEDs, or any spectral region that contains a significant contribution from predicted lines. 
However, progress can be made by making corrections to these high resolutio n libraries for 



the ef fect of predicted lines on the broadband SED, as for example done in ICoelho et al 
( 20071 ). 



How Accurate are the Theoretical Libraries? Normally, observed stars are assigned 
physical parameters based on comparison to models. However, if one wants to constrain 
the models with observed stellar data then one requires an independent estimate of the 
physical properties of stars. This is a significant obstacle to assessing the accuracy of 
the models. One can obtain essentially model-independent estimates of T e g via angular 
diameter meas urements, but the s e are very difficult to obtain for giants, let alone main 
sequence stars (jPerrin et al.lll998l . iBovaiian et al.ll2012a). A more widely used technique is 
the infrared flux method (IRFM), first proposed bv lBlackwell fc Shallisl (jl977V ). The basic 
idea is that the flux of a star depends only weakly on T e ff in the IR, and so one can deduce 
angular diameters from fluxes alone with only a weak model dependence. This technique 



underpins nearly all existing co l or-T e ff relations (e.g., . 
19961 ; I Ramirez fc Melendezj [20051 ; ICasagrande et al.ll201oT ) 



Alonso. Arribas fc Martinez-Roger 



The latest generation of m odel spectra succeed in reproducing the broadband colors for 
FGK dwarfs and warm gian ts ( Bertone et al.|[20o3 . Kucinskas et al.ll2005l . Martins fc Coelhol 
2007llBovaiian et al.lE)12bft . Not surprisingly, the models continue to have difficultly fitting 
the SEDs of the coolest stars and the flux in the Wien tail of the flux distribution (typically 
< 5000 A for types G and later) . The latter difficulty arises because the Wien tail is 
extremely sensitive to Teff and so even minor changes to the model assumptions lead to large 
changes in that spectral region. In addition, metal line blanketing becomes very strong in the 
blue/UV, and so the requirements on the line lists are demanding. IMartins fc Coelhol ( 20071 ) 
explored the ability of three modern synthetic libraries to reproduce a variety of spectral 
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Figure 2: 

Physical properties of the MILES empirical spectral library, separated by [Fe/H] (solid symbols). 
Isochrones from the Padova group are overplotted for i = 3x 10 6 — 10 1 yr. Typical errors on T e fj 
and log g are 60 — 100 K and 0.2 dex, respectively iCenarro et al.ll2007h ■ Notice that the MILES 
library covers the lower main sequence and RGB both at high and low metallicity, but it only 
sparsely covers the upper main sequence and supergiants. This figure highlights the difficulty in 
constructing SPS models based on empirical stellar libraries, especially at low metallicity. 



features of observed stars. Overall the agreement was found to be satisfactory, although 
there were notable deficiencies. Some models were unable to reproduce the observed CH 
features in cool stars, while other models failed to match the observed Mgfc and MgH 
features. The majority of the model shortcomings were restricted to stars with T c g < 4500 
K. M odern models also have difficulty reproducing the streng th of the TiO bands in cool 
stars ( Kucinskas et al. 2005 , Allard. Homeier fc Frevtaa 2011 ) . 

The situation regarding the H2O band strengths in the NIR should serve as a caution- 
ary tale: models have for decades predicted water bands in cool stars that are stronger 
than observations. Part of the problem seems to have been the mo l ecular line lists, which 
improved steadily over the years, culmin ating in the iBarber et all (200fj|) line list, which 
contains 500 million transitions. However, Allard. Homeier 



^revtag 



( 2011 ) demonstrated 

that the recent revision of the solar abundance scale bv lAsplund et al.l (|2009T ) . which entails 
a factor of two reduction in the oxygen abundance, has a much larger effect on the predicted 
strength of the water bands. Their new models with updated solar abundances reproduce 
remarkably well the NIR colors of M dwarfs, which have strong water band features. Thus, 
not only the atomic and molecular parameters but also the abundance patterns can have a 
very significant effect on theoretical spectra. 



Empirical Libraries The strengths and weaknesses of the theoretical libraries are the 
weaknesses and strengths of the empirical libraries. Empirical spectra of course do not suffer 
from issues with line lists, treatment of convection, etc., but they are plagued by standard 
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observational constraints such as correction for atmospheric absorption, flux calibration, 
and limited wavelength coverage and spectral resolution. Worse, the empirical libraries are 
woefully incomplete in their coverage of parameter space. This is a long-standing issue that 
is difficult to address owing to the fact that empirical libraries are drawn from samples of 
stars in the solar neighborhood. For example, hot main sequence stars at low metallicity 
are very rare, as are stars in rapid phases of stellar evolution such as WR and TP-AGB 
stars. 

One of the first com prehensive empirical stellar spectral libraries was constructed by 
Gunn fc Strvkej dl983T). Later optic al/NIR libr aries included I Pickles! (|l998l ). Jo nes (1999), 
ELO DIE dPrugniel fc SoubirarfcooJ). STELIB jLe Borgne et al.ll2003h. Indo-US Jvaldes et al 



2004 ). NGS L Jcregg et al.ll2006l. iHeap fc Lindleril201ll l. MILES JSanchez-Blazauez et al 
2QM), IRTF l|Ravner. Cushing fc Vaccal2009h . and the X-shooter library (jchen et all2011n" 



Othe r, specialized libraries include an ultraviolet atlas compiled from IUE data ijFanelli et al 



stellar library (|Ardila et all |201cF . The Sloan Digital Sky Survey (SDSS) has obtained 



19921 ). a TP-A GB library jLancon fc Mouhcinelboool ). and the Spitzer Space Telescope IR 



spectra of several hundred thousand sta rs and has devo ted considerable effort to measur- 
ing stellar parameters from the spectra ( Lee et al. 20081 ) . As yet there is no official stellar 
spectral library based on the SDSS stars, but it would clearly provide a great leap forward 
both in terms of coverage in parameter space and uniformity in spectral quality. 

An example of the difficulties posed by empirical libraries for SPS model construction 
is shown in Figure [2] This figure shows the location in the HR dia gram of all sta r s from 
the MILES spectral library, with stellar parameters determined bv lCenarro et al.l (|2007l ). 
MILES is the empirical library today that provides the greatest coverage in terms of logg, 
T e ff, and [Fe/H] with a total of 985 stars. In this figure the metal-rich and metal-poor 
stars are split into two panels. Overplotted are isochrones from the Padova group for ages 
ranging from 10 6 ' 6 — 10 1 yr. While the lower main sequence and red giant branch are 
well covered by the empirical library over a wide range in metallicity, hotter stars are rare, 
especially at lower metallicity. Because the HR diagram is covered so sparsely in this regime, 
constructing SPS models based on empirical libraries at young ages (< 1 Gyr) is clearly 
challenging. 

A second, related problem with empirical libraries is irregular coverage in the HR di- 
agram. Interpolation w i thin t he library, necessary for SPS model construction, can be 
difficult. IVazdekis et~aH l|20ld ) has developed a sophisticated algorithm to deal with this 
issue. These authors assign quality numbers to the resulting SPS models based on the num- 
ber of empirical spectra comprising each model. This is a profitable approach that should 
be standard in the field because it allows users to assess the reliability of various regions of 
parameter space. 

A third problem associated with empirical libraries is the assignment of physical param- 
eters to stars, including log<?, T e g, and [Fe/H]. Of particular note is the T c a determina- 
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Percival et al 



tion, which carries an uncertainty of order 100 K ( Alonso. Arribas 
Ramirez fc Melendezl20o"il ; ICasagrande et al.l201O . Using theoretical models 
( 20091 ) explored the impact of uncertainties in T e g, logg, and [Fe/H] for individual stars on 
SPS models. They found that for models with old ages (> 4 Gyr), the effect of changing the 
temperature scale by 100 K had significant effects (0.1 — 0.4 Ain EW in the most extreme 
cases) on a variety of spectral absorption features, including the hydrogen Balmer lines and 
several commonly employed iron and magnesium lines used to interpret the stellar popula- 
tions of early-type galaxies. The work of Percival & Salaris was mostly exploratory, and it 
would be desirable to consider the full propagation of these sorts of uncertainties into the 
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Figure 3: 



Variation in the spectrum of a 13 Gyr SSP due to changes in individual elemental abundances. All 
abundance changes are +0.3 dex except for C and Ca which are varied by +0.15 dex. Spectra 
have been broadened to a velocity dispersion of 150 kms -1 . Notice the different x and y-ax.es in 
each panel. Prominent spectral features include the CH, CN, NH, C2, Ca I, Ca II H&K, Ca II 
triplet, MgH, and Mg b (Mg 1) features, along with numerous atomic Fe and Ti features. 
Computed from the models of lConrov fe van Do kkum ( 20 12al ). 



derived properties of galaxies. 

Yet another issue with empirical libraries is the abundance patterns of the stars. It is 
well-known that low metallicity stars in the Galaxy tend to be a— enhanced, such tha t 
[Mg/Fe]« 0.0 at [Fe/H]« 0.0 and [Mg/Fe] « 0.4 at [Fe/H]< -1.0 l|Edvardsson et al.lll99ot ). 



Thus, any model based on empirical stars at low metallicity must somehow correct for th is 
[a/Fe] 'bias' in the models (see e.g.. iThomas. Maraston fe Bendenl2003bl ; ISchiavonll2007l ). 



Variable Abundance Patterns Moderate resolution spectra (i.e., R ~ 1000 — 5000) 
contain a wealth of information on the detailed abundance patterns of the stellar popu- 
lations. In the context of SPS models, abundance ratios have historically been estimated 
through the analysis of spectral indices. Ideally, an index is sensitive to a single feature 
(such as an atomic absorption line or a molecular band head). One usually defines a fea- 
ture bandpass and one or two pseudocontinua, and then an equivalent width (EW), or 
in some cases a magnitude, can be measured. The Lick/IDS system is th e most popular 



index system, defining 21 indices in the wavelength range 4000 — 6500 A i Burstein et al 



1984 



(e.g.JFanelli et al. 



Worthev et al.lll994h. Other index systems have been defined at other wavelengths 
19921 . 1 Alvarez et al.ll200a ICenarro et "ZlboOll . IServen. Worthev fe Briiev 



2005). In practice, the use of indices is greatly complicated by the fact that there are rarely 
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if ever clean regions of the spectrum from which one can estimate the continuum level. The 
strength of each index thus depends not only on the feature of interest but also the features 
in the pseudocont i nua. 



Tripicco fe Belli (|1995T ) were the first to assess the sensitivity of the Lic k/IDS indices to 



variat ion in the abundances of 10 elements with theoretical spectral models. iKorn. Maraston fc Thomas 
( 2005 ) provided an update to these 'response f unctions' with updated line l ists an d transition 



probabilities, and for a range in metallicity. IServen. Worthev fc Brilevl (|2005l ) considered 



variation in 23 elements, including several neutron-capture elements. These authors defined 
new indices sensitive to a host of elements not considered in the earlier work. The Serven 
models were based on theoretical spectra of only two stars per abundance pattern ; a turnoff 



star an d a luminous giant, and the model effectively had a fixed age of ~ 5 Gyr. iLee et al 



( 2009al ) provide d response functions for Lick in dices by employing a larger number of syn- 



thetic spectra. IConrov fc van Dokkuml |2012ah produced theoretical stellar spectra with 



variation in 11 elements. These models contained theoretical spectra at 20 points along 
the isochrone, for ages from 3 — 13.5 Gyr. Rather than focusing on spectral indices, these 
models made predictions for the response of the full spectrum from 3500 A— 2.4 /im to abun- 
dance variations. Sample response spectra from that work are shown in Figure In this 
figure the relative response of the spectrum to an increase in the abundance of a particular 
element is shown. The full spectrum is rich in diagnostic features beyond what any index 
system can capture. 

To construct models with arbitrary abundance patterns one would want to create syn- 
thetic spectra for each possible abundance pattern. When considering large numbers of 
elements, such an approach becomes computationally infeasible. Instead, the standard tech- 
nique is to create arbitrary abundance patterns by treating t he effect of each element on the 
spectrum as being independent of the other elements (e.g., Thomas. Maraston fc Bender 
2003bl ; ISchiavorj2007l : ILee et alj|2009al ; IConrov fc van Dokkurrj2012ah . This is a reasonable 
assumption for trace elements, but it is less realistic for elements such as C, N, O, Na, Ti, 
and Fe, which affect the opacity, free electron density, and/or molecular equilibrium. There 
has been surprisingly little work in the literature testing this assumption. 



2.1.4 The IMF The initial distribution of stellar masses along the main sequence, 
know n as the stellar initial ma s s function ("IMF) , has been studied extensively for decades 



known as trie stellar initial mass function (1M1), lias been studied extensively lor decade; 
(e.g. . ISalpeteill 19551 . [jcaloll 198rj. IScalo et al.lll998l. lKroupal|200ll , Ichabrierll2003n . As empha- 



sized in a recent review bv lBastian. Covey fc Meve il l|2010h . there is no compelling evidence 
for variation in the IMF from direct probes, e.g., star counts. The canonical Salpeter IMF 
has the form dN/dM oc M~ x with x = 2.35. The IMF measured in the solar neigh borhood 
deviates from the Sa lpeter form only at M < 1 Mq , where x becomes shallower (jKroupa 
200ll . IChabrierlliool ) . 



From the perspective of SPS, the IMF (1) determines the overall normalization of the stel- 
lar mass-to- light ratio, M/L; (2) determines the rate of luminosity evolution for a passively 
evolving population; (3) affects the SED of composite stellar populations; (4) has a small 
effect on the shape of the SED of single stellar populations. The last point is due to the fact 
that the integrated light of a coeval population is overwhelmingly dominated by stars at 
approximately the same mass, i.e., the turnoff mass. (4) will not hold for IMFs that depart 
dramatically from the Salpeter IMF. (3) arises because composite populations have a range 
of turnoff masses that contribute to the integrated light. The fractional contribution of 
stars of various masses to the total number of stars, stellar mass, and bolometric luminosity 
is shown in Figure [4] This figure demonstrates quantitatively that low mass stars dominate 
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Figure 4: 

Fractional contribution to the total number, mass, and bolomctric luminosity as a function of 
stellar mass for a 13.5 Gyr solar metallicity model. Lines correspond to different IMFs: a 
bottom-heavy with logarithmic slope x = 3.0 (blue line); Salpeter (x = 2.35; black line); MW IMF 
(specifically a Kroupa IMF; red line); a bottom-light IMF (specifically of the form advocated by 
van Dokkum (2008); green line). The inset in the right panel shows the cumulative luminosity 
fraction in logarithmic units. Low mass stars dominate the total number and mass in stars, but 
contribute a tiny fraction of the luminosity of old stellar populations. 



the stellar mass and number of stars in a galaxy, but contribute only a few percent to the 
bolometric light of an old stellar population. At younger ages the light contribution from 
lo w mass s tars is even less. 

Tinslevl (|l980h demonstrated that the evolution in M/L for a passively evolving stellar 



population is sensitive to the logarithmic slope of the IMF, x, at the main sequence turnoff 
point. The logarithmic evolution of the luminosity per logarithmic time (dlnL/dlnt) scales 



Dokkuml200i . lConrov, Gunn fc White 



linea rly with x, at least for plausible values of x (see also|; 
2009). This dependency arises because the giant branch dominates the luminosity for all 



plausible values of x, and so the IMF determines the rate at which the giant branch is 
fed by turnoff stars. Steeper IMFs imply that the giant branch is more richly populated 
with time, and therefore the natural luminosity dimming is reduced. For sufficiently steep 
IMFs (e.g., x > 5), the unevolving dwarfs would dominate the light, and so the integrated 
luminosity would be approximately constant over time. 

It is somewhat less- well appreciated that the IMF above 1 Mq also strongly affects the 
shape of the SED of composite stellar populations. In composite populations the SED is 
influenced by stars with a range of masses (see Section 1 2. 3[) . and so the IMF must play an 



important role (see e.g.. lPforr. Maraston fe Toninill201 



pi), 



2.2 Dust 

Interstellar dust is a component of nearly all galaxies, especially those that are actively 
star- forming. Dust plays a dual role in SPS, both as an obscurer of light in the UV through 
NIR and as an emitter of light in the IR. For both historical and theoretical reasons these 
two aspects are often modeled independently of one another. Indeed, these two aspects 
are sensitive to different properties of a galaxy (e.g., dust obscuration is highly sensitive to 
geometry while dust emission is more sensitive to the interstellar radiation field), and so it 
is not unreasonable to decouple the modeling of these two components. 
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2.2.1 Attenuation Dust grains obscure light by both absorbing and scattering starlight 
From observations of individual stars one can infer the total extinction along the line of sight 
by comparing an observed spectrum to the expected spectrum of the star in the absence of 
dust (the latter is typically obtained from models in conjunction with an estimated T c s and 
log g of the source). Thus, extinction measures the total loss of light along a single line of 
sight. Observations in the Milky Way (MW) and Large and Small Magellanic Clouds (LMC 
and SMC, respectively) have been obtained for a number of sightlines, enabling construction 
of average extinction cur ves for these galaxies ( Cardelli. Clayton fc Mathisl 19891 ; Pei 1992 ; 
Gordon fc Clayton 1998h . In the MW and LMC the only striking feature in the extinction 
curve is the broad absorption at 2175 A. This feature is absent in three out of four sightlines 
to the SMC, and it is weaker in the LMC than in the MW. The grain population respon- 
sible for this feat ure is not kn own, but polycyclic aromatic hydrocarbons (PAHs) are a 
leading candidate ( Drainei 2003 ) . In general, dust extinction is a co nsequence of the optical 
prop erties of the grains and the grain size and shape distribution ( Weingartner fc Drainei 



2001 



When modeling SEDs of galaxies, the relevant concept is dust attenuation, which differs 
from extinction in two important respects: (1) light can be scattered both out of and into a 
given line of sight; (2) the ge ometrical distr ibution of dust with respect to the stars strongly 
affects the resulting SED (see CalzettilbOQl . for an extensive discussion of these issues). The 
total dust attenuation in a galaxy can be estimated by analogy with how one estimates dust 
extinction: a spectrum of a galaxy is obtained and compared with the expected spectrum 
of the same galaxy in the absence of dust. For obvious reasons, estimating dust attenuation 
is considerably more complex than estimating dust extinction. 

Although the shape of the dust attenuation curve depends on the star-dust geometry, 
grain size distri bution, etc., in a com plex manner, several general rules of thumb can be 
stated (see e.g., Witt fc Gordon 2000l . for a more thorough discussion). The simplest dust 
geometry, that of a homogenous foreground screen, yields an attenuation curve whose shape 
depends only weakly on the total dust column density (the weak dependence arises from the 
scattered light component). More complex geometries generally yield attenuation curves 
that become greyer (i.e., shallower) as the column density increases. Clumpy interstellar 
media also result in greyer attenuation curves than their homogenous counterparts. In all 
cases, the total attenuation optical depth is less than the amount of extinction that would 
be produced by the same column density because some scattered light is always returned 
to the line-of-sight. Finally, because the 2175 A dust feature is believed to be due to pure 
absorption, the effects of radiative transfer will cause this feature to respond differently to 
geometrical effects compared to the rest of the attenuation curve. The use of a single dust 
attenuation curve for analyzing a wide range of SED types is therefore without theoretical 
justification. 

In practice, most SPS modelers include dust attenuation by fixing the shape of the at- 
tenuation curve and fitting for the normalization. Popular attenuation curves include the 
Calzetti law, the MW, LMC, and SMC extinction curves, or the time-dependent attenu- 
ation model from I Chariot fc Falll (|2000l ). The impact of the chosen attenuation curve on 
derived properties of galaxies will be discussed in later sections. 



2.2.2 Emission Beyon d A 
emission from dust grains 



10 jim the SEDs of normal galaxies are dominated by 
Mathis. Rumpl fc Nordsieckl (1977) were the first to postulate 
that the dust grain population is comprised of both silicate and carbonaceous grains. In 
modern theories, the carbonaceous grains are assumed to be PAHs when they are small, 
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and graphite when they are large (jDraind 120031 ). The observed dust emission spectrum 
results from exposing these grains to a range of interstellar radiation field strengths. 

A variety of models exist that combine grain size distributions and grain optical prop- 
erties with models for starlight ( or simply the radiation field) to predict IR emission 
dDesert, Boulanger fc Pugetlll990l; ISilva et aljfliii; bevriendt. Guiderdoni fc Sadatlll99E 



Popescu et al.ll2000l ; Dale et al.ll200ll;|Piovan, Tantalo fc Chiosill2006l ; |jonssonl2006l ; lDraine fc~ 



20071 ; iGroves et al.l I 2008 ; Popescu et al.ll201lh . The modeling of PAH emission features 
(the most prominent being at 3.6 fim, 6.2 fj,m, 7.7 jim, 8.6 Atm ; and 11.3 /xm) has become 



drama tically more sophisticated since the early attempts by Desert. Boulanger fc Puget 
( 1990l ). culminating in state-of-the-art models bv lDraine fc Lil (|2007T ) At long wavelengths 
(A > 50 /jm) the emission is dominated by grains at a nearly steady temperature of ~ 15 — 20 
K and contributes ~ 2/3 of the total IR luminosity. At shorter wavelengths the IR emis- 
sion arises from single photon heating of dust grains (including PAHs) and accounts for 
the remaining ~ 1/3 of the total IR emission (see the review bv lDraine|[2003l . for details). 
The IR emission at the shortest wavelengths (A < 12 (im) is supplied almost entirely by 
PAHs. In detail these relative contributions will depend on the grain composition and size 
distribution and the interstellar radiation field. 

The models listed above are not always well-suited for interpreting large numbers of ob- 
served IR SEDs because they contain a large number of parameters, require knowledge of 
the star-dust geometry, and/or require radiative transfer calculations. Simpler models for 
dust emissio n have therefore been d evelo ped in para l lel to the more complex ones. The IR 



templates of lcharv fc Elbazl |200ll ) and bale etaH |200ll ) are widely used for estimating 
bolometric luminosities and fc— corrections. These templates are based on sophisticated 
models, but are constrained to match observations of normal star-forming and starburst 
(i.e., ULI RG) galaxies. The resu lting templates are functions of only one variable. For ex- 
ample, the Chary fc Elbazl ( 2001 ) templates are a function of the bolometric IR luminosity. 



da Cunha. Chariot fc Elbaz ( 20081 ) have developed a simple phenomenological model for 
dust emission that consists of a series of modified blackbodies for the thermal dust emission 
and for the emission from stochastically heated dust grains. In addition, they include an 
empirical spectrum of M17 to represent the PAH emission. 



2.2.3 Dust Around AGB Stars Stars experience very high mass-loss rates as they 
climb up the AGB (as high as 10~ 4 Mp) yr - 1 in the superwind phase). The mass lost is 



often observed to be dust rich (|Bediinlll987l ). These stars are often heavily dust obscured 
in the optica l and emit co piously in the IR; they are often so dusty that their SEDs peak 
at ~ 10 fim (|BediirJll987l ). Dust around AGB stars is important from the standpoint of 
SPS for two reasons: it will diminish the importance of AGB light in the optical and NIR, 
and it will contribute addition al flux in the mid-IR be yond what would be expected from 
standard dust models (see e.g., Kelson fc Holden 20101 ) . 



Despite the obvious importance of AGB d ust, this aspect is included in fe w SPS mod 



els. N o table exceptions include the mod els of |Bres^aii J _Granato & Silva 



Silva et al 



( 19981 ). IPiovan. Tantalo fc Chip's] (|2003l ). and iGonzalez-Lopezlira et alj (|2010l ) , which are 



all based on theoretical models for AGB stars, their mass-loss rates, and dust formation and 
composition. Even the use of currently available empirical libr aries of AGB stellar spectr a 
require reddening the spectra to account for circumstellar dust ( Lancon fc Mouhcine 20021 ) . 
Increasingly sophis ticated radiative transfer models of dusty circumstellar envelopes are be- 
ing d eveloped (e.g., Groenewegerj|2o"l2l Sargent. Srinivasan fc Meixnerl201ll . lSrinivasan. Sargent fc Meixner 
201ll ). The outputs from these models should be incorporated into SPS codes as a standard 
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Figure 5: 

Top Panels: Fractional contribution to the total flux from stars in various evolutionary phases, for 
three different SFHs. The left panel is representative of a galaxy that formed nearly all of its stars 
very rapidly at early times, the middle panel is representative of a typical star-forming galaxy at 
z ~ 0, and the right panel may be representative of the typical galaxy at high redshift. Flux 
contributions are at 13 Gyr (solid lines) and 1 Gyr (dashed lines) after the commencement of star 
formation; all models are sola r mctallicity, dust-free, and are from FSPS (v2.3; 
IConrov, Gunn & White 2009). Labeled phases include the main sequence (MS), red giant branch 
(RGB), asymptotic giant branch (AGB, including the TP-AGB), post-AGB (pAGB), and the blue 
and red horizontal branch (bHB and rHB). Bottom Left Panel: Fractional flux contributions for 
stars more massive than 20 Mq and 60 Mq for the SFH in the middle panel of the top row. 
Bottom Right Panel: Light-weighted age as a function of wavelength for the same SFH. The 
dashed line indicates the corresponding mass-weighted age. 

practice. 

2.3 Composite Stellar Populations 

The simple stellar populations discussed in Section 12.11 are the building blocks for more 
complex stellar systems. Composite stellar populations (CSPs) differ from simple ones in 
three respects: (1) they contain stars with a range of ages given by their SFH; (2) they 
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contain stars with a range in metallicities as given by their time-dependent metallicity 
distribution function, P(Z,t); (3) they contain dust. These components are combined in 
the following way: 



/csp(t) = / / \SFR(t~t')P{Z,t-t')fssp(t',Z)e- Td< - t ' ) +Af duBt {t',Z) J dt' dZ 
Jt'=o J z=o \ ) 

(2) 

where the integration variables are the stellar population age, t' , and metallicity, Z. Time- 
dependent dust attenuation is modeled via the dust optical depth, Td(t') and dust emission 
is incorporated in the parameter /dust- The normalization constant A is set by balancing 
the luminosity absorbed by dust with the total luminosity re-radiated by dust. 

The SFH can in principle be arbitrarily complex, although simple forms are usually 
adopted. By far the most popular is the exponential, or r— model, where SFR cx e~ t//r . This 
form arises natur ally in scenario s where the SFR depends linearly on the gas density in a 
closed-box model ( Schmidtll95S ). Recently, rising SFHs have become popu lar to explain the 
SEDs of high-redshift galaxies ( Maraston et al. 20ip| . Papovich et al ] l201lh . Rising SFHs at 
early times seem to be a natural consequence of galaxy evolution in a hierarchical universe 
( Finlator. Dave fc Qppenheimerll2007l : iLee et al.ll2010f ). Functional forms that incorporate 
an early phase of rising SFRs with late-time decay, such as SFR oc e~ t//r , may therefore 
become more popular amongst modelers. 

The treatment of metallicity in composite stellar populations is usually even more simplis- 
tic than the treatment of SFRs. The widely adopted simplification is to replace P(Z, t) in 
Equation [5] with a 8— function. In other words, a single metallicity is assumed for the entire 
composite population. The impact of this simplification on the SPS modeling procedure 
has not been extensively explored. 

The approach to Equation [2] outlined above is standard but not universal. A notable 
exception is the technique of fitting non-par ametric SFHs and metallici t y histories in either 



a pre-defined or a daptive set of age bins ( Cid Fernandes et al. 2005 . Ocvirk et al. 20061 . 



Toieiro et ail 120091 ). The latter approach is computationally expensive since as many as 
thirty parameters are simultaneously fit to the data. Very high quality data are also a 
prerequisite for non-parametric techniques. Nonetheless, they offer the promise of less 
biased reconstruction of the SFH and metallicity history of galaxies based solely on their 
SEDs. 

Figure [S] presents some basic properties of CSPs. The top panels show the fractional flux 
contribution from stars in different evolutionary phases for three representative SFHs. At 
late times the RGB and red HB dominate the red and NIR flux of r-model SFHs, as is well- 
known. However, at young ages and/or for rising SFHs, the red and NIR are also influenced 
strongly by AGB stars, suggesting that the NIR can be susceptible to large uncertainties 
(due to uncertainties in the modeling of the AGB, as discussed in detail in later sections). 
The lower left panel shows the fractional flux contribution from stars more massive than 
20 Mq and 60 Mq for a r = 10 Gyr SFH. Massive star evolution is uncertain due to the 
complicating effects of rotation and binarity, and so this panel provides a rough sense of the 
extent to which uncertainties in massive star evolution will affect SED modeling. Finally, 
the lower right panel shows the light-weighted age as a function of wavelength, again for 
a t = 10 Gyr SFH. The dashed line indicates the mass-weighted mean age. Already by 
0.5 )im (approximately the V— band), the light- weighted age reaches its maximal value, 
again suggesting that going to redder restframe wavelengths is not providing substantial 
new information on the SFH. Notice also that the maximum light-weighted age never reaches 



IS 



the mass- weighted age. We return to this point in Section [4] 



2.4 Nebular Emission 



Although the effects of nebular emission on SEDs will not be discussed in detail in this 
review, a brief overview of this component is given for completeness. 

Nebular emission is comprised of two components: continuum emission consisting of 
free-free, free-bound, and two photon emission, and recombination line emission. Several 
photoionization codes exist that make predictions for the nebular emi ssion as a function of 
the physical state of the gas, includi ng CLOUDY (jFerland et alJll99Sl ) and MAPPINGSIII 
( Groves. Dopita fe Sutherland! hood ). Other approaches can be taken t o the modeling of 



nebular emission lines. For example, lAnders fc Fritze-v. Alvensleben (2003) implement 



non-hydrogen emission lines based on observed line ratios as a function of metallicity. The 
nebular emission model must then be self-consistently coupled to a model for the starlight. 
Several groups have done this, with varying degrees of complexity, and with some only 
including nebu l ar continuum, others only line emission, and others including both (e.g. , 



Leitherer et alj|l999l. ICharlot fc LonghettilbOQll, IPanuzzo et alj|2003l . I Groves et Zll2008l . 



Molla. Garcia- Vargas fe Bressan 20091 . Schaerer fe de Barrosl 201C ). 

The effect of nebular emission on the SED is complex, especially when line emission is 
considered. As a general rule, nebular emission is more important at low metallicity and 
at young ages. In such cas es the contribution of nebular emissi on to broadband fluxes 
can be as high as 20 — 60% ( Anders fe Fritze-v. Alvensleben! 2003 ). Nebular emission will 
also be more important at high redshift because a feature with a fixed restframe EW will 
occupy a larger fraction of the filter bandpass due to the redshifting of the spectrum. The 
effect of nebular emission therefor e cannot be ignored at high redshift, where hi gh SFR, low 
metallicity galaxies are common (|Schaerer fe de Barros!l2010l.lAtek et al.ll201lh. 



It is noteworthy that amongst the most widely used SPS codes, incl uding lBruzual fe~ larlo^ 



Marastor: 



tiy tnat amongst tnc most widely used al^a codes, including cruzual &z o larlot 
3 J2005I). PEGASE llFioc fe . Rocca-Volmerangell99"7h , STARBURST99 l|Leitherer et al 



19991 ). and FSPS (jConrov, Gunn fe White! 120091 ). only PEGASE and STARBURST99 



elude nebular continuum emission, and only PEGASE includes both nebular continuum and 
line emission. Since nebular emission is relatively straightforward to implement (notwith- 
standing assumptions about the physical state of the gas), it should be a standard compo- 
nent of all SPS codes. 



2.5 Fitting Models to Data 

The SPS models described in this section are most frequently used to measure physical 
parameters of stellar populations. This is achieved by fitting the models to data, either in 
the form of broadband SEDs, moderate-resolution optical/NIR spectra, or spectral indices. 
The SSPs are usually taken as given, and the user then fits for a variety of parameters 
including metallicity, dust attenuation, and one or more parameters for the SFH. If IR 
data are available, then additional variables must be considered. The fitting techniques 
vary but are generally limited to grid-based x 2 minimization techniques. As the number 
of para meters increases Markov Chain Monte Carlo techniques become in creasingly more 



efficient ( Conrov. Gunn fc Whitel2009l ; Acquaviva. Gawiser fc GuaitalboTl ). One must also 



be aware of the influence of the chosen priors on the derived parameters; in cases where 
parameters are poorly constrained the prior can have a significant effect on the best-fit 
values (see e.g.. lKauffmann et al . 2003, ISalim et al.ll2007l . lTavlor et al.ll201ll . for discussion). 
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Moreover, it is highly advisable to derive the 'best-fit' parameters from the marginalized 
posterior distribution, rather than from the minimum of y 2 , since the likelihood surfa ce can 
often be highly irregular (see e.g., Bundv. Ellis fc Conselicg 20051 ; Taylor et al. 2011). 

When fitting SEDs it is important to remember that only the shape is being used to con- 
strain the model parameters. In other words, the SED shape (or, more crudely, broadband 
colors) constrains parameters such as M/L, specific SFR (SSFR=SFR/M), dust attenua- 
tion and metallicity. To obtain stellar masses one needs to multiply M/L by the observed 
luminosity and to obtain SFRs one then multiplies SSFR by M. This also holds true 
for M/L ratios and SSFRs inferred from EWs of emission lines, spectral indices, and full 
spectral fitting. 

The reader is referred to IWalcher et alJ (|201ll ) for further details on SED fitting tech- 
niques. 



3 MASS-TO-LIGHT RATIOS & STELLAR MASSES 
3.1 Techniques &: Uncertainties 

There are three basic techniques for estimating the stellar mass-to-light ratio, M/L, of a 
galaxy: (1) using tabulated relations between color and M/L; (2) modeling broadband 
photometry; (3) modeling moderate resolution spectra. The first technique is the simplest 
to use as it requires photometry in only two bands and no explicit modeling. The other 
techniques require construction of a library of models and a means to fit those models to 
the data. How do these different techniques compare? 

3.1.1 Color-based M/L Ratios The c olor-based M/L estimators have their origin in 
the pioneering work of Bell fc de Jonj (2001). These authors used a preliminary version of 



the lBruzual fc Charloti |2003l ) SPS model to chart out the relation between M/L and color 



as a function of metallicity and SFH. Remarkably, they found that variation in metallicity 
and SFH (parameterized as a r model) moved galaxies along a well-defined locus in the 
space of M/Lb vs. B — R, suggesting that the B — R color could be a useful proxy 
for M/ Lb- Perhaps most importantly, they demonstrated that the dust reddening vector 
was approximately parallel to the inferred color-M/L relation, implying that dust should 
have only a second-order effect on derived M/L values. They also investigated NIR colors 
and found much larger variation in M/L at fixed I — K, due mostly to variation in the 
SFH. Late bursts of SF complicate the interpretation of color- M/L relations by driving the 
M/L ratios lower at fixed color compared to smoothly declining SFHs. Their analysis was 
for a fixed IMF; allowing for IMF variation will change the M/L while leaving the color 
basically unchanged. These authors concluded that for a fixed IMF and assuming that 
large starbursts and low-metallicity systems are not common properties of galaxies, and 
neglecting uncertainties in the SPS models, one could estimate M/L from a single color to 

an accuracy of ~ 0.1 — 0.2 dex. 

In subsequent work. lBell et al. I (|2003h analyzed the optical-NIR SEDs of ~ 12, 000 galaxies 



with SDSS and 2MASS photometry. They constructed a grid of SPS models with varying 
SFHs, including bursts, and metallicity. They did not allow for reddening due to dust. They 
derived best-fit M/L values from SED fitting and created new, observation ally-constrained 



color- M/L relations. The resulting B — Rvs. M/Lb relation was similar to lBell fc de Jone 



(2001), but the NIR relation was considerably more shallow than the earlier work. At the 



bluest colors, the M/Lk ratios differed by ~ 0.3 dex between the old and new relations. The 
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difference arose due to a popul a tion o f galaxies with blue observed colors and high inferred 



M/L ratios, which iBell et all |2003h interpreted as being due to the lower metallicities 



allowed in their fitting procedure compared to th e analytic model s of Bell & de Jong. 
These differences are of some concern because the Bell et al. ( 20031 ) color-M/L relations 
a re very widely used to d e rive ' cheap' stellar mass estimates of galaxies. 

Zibetti. Chariot fc Rix (|2009h revisited this issue with new SPS models from Chariot 



Bruzu al (in prep). They followed a somewhat different approach compared to IBell et al 



( 20031 ) in that they constructed color-M/ L relations directly from a library of model galax- 
ies. This approach is sensitive to how one chooses to populate parameter space, since 
the resulting col or-M/L relations are simple averages over the model space. In contrast, 



Bell et al. 1 2003h constructed color-M/L relations only from those models that provided 



good fits to observed optical-NIR photometry. Zibetti et al. allowed for dust attenuation, 
and they also allowed for larger mass fra ctions of stars formed in recent bursts compared to 
Bell et al. Izibetti. Chariot fc Rixl |2009l) found substantial differences in the g — i color-M/L 
relation compared to Bell et al. (~ 0.5 dex at the bluest colors), which they attributed 



to the ir allowance for much s tronger bursts of star formation than allowed in IBell et al 



( 20031). iGallazzi fc Belli (|2009h also found a very different color-M/L relation compared 



Bell et al 



to 



J2003I). wh ich th ey attributed to the diffe rent age di s tribu tions in their model 



library. iTavlor et al.l (|201ll ) followed the philosophy of lBell et al. I (|2003l ) in fitting models to 
observed SEDs in order to derive color-M/L relations. They found systematic differences 
of order 0.1 — 0.2 dex compared to the Zibetti et al. relations, which they attributed mostly 
to their decision to select only those models that fit the photometric data, and also the 
tr eatment of dust redden i ng. 

Zibetti. Chariot fc Rixl (|2009l ) found a large difference in M/Lh ratios esti mated from i — 

H colors between the new Chariot & Bruzual models and the previous version (jBruzual fc Chariot 



20031 ). with typical offsets of 0.2 — 0.4 dex. This difference is largely due to the different 
treatment of the TP-AGB evolutionary phase (see Section 13.21 for details). In contrast, 
optical color- M/L relations showed li ttle difference be tween the two SPS models, and also 
showed slightly less discrepancy with Bell et al. 1 20031 ). These authors argued that optical 
color-M/L relations should generally be preferred to NIR-based ones owing to the larger 
impact of the uncertain TP-AGB phase in the NIR. Exceptions to this rule could be made 
for galaxies with very large amounts of dust reddening and/or very strong starbursts, in 
which cases the optical SED becomes a poor indicator of the SFH and hence M/L. The up- 
shot is that NIR color-based M/L estimates should be treated with caution, since they are 
subject to large stellar evolution-based systematics. However, both optical and NIR-based 
M I L estimates will depend on the allowed range of SFHs and metallicities in the model 
library, with large systematics appearing for the bluest colors. 

3.1.2 M/L From Broadband and Spectra l Fitting Technique s The modern era 
of SED fitting with SPS models was ushered in by Sawicki fc Yee ( 19981 ) and Giallongo et al 



(1998), who a nalyzed optical-NIR b roadband photometry of high-redshift galaxies. The 
model grid of Sawicki fc Yeel ( 19981 ) consisted of a r— model SFH with varying r and 
start time for the SFH, and variation in metallicity and redd ening. Subsequent model- 



ing of broadband data has largely followed this approach (e.g., Brinchmann fc Ellis! 2000l . 



Papovich. Dickinson fc Ferguson 2001 . Shaplev et al. 2001 . Salim et al. 20071 ). Genericallv. 



when fitting broadband SEDs one finds that stellar masses for 'normal' galaxies (i.e., not 
including pathological SFHs) can be recovered at the ~ 0.3 dex level (l<r uncertainty). 
This uncertainty does not include potential systematics in the underlying SPS models; 
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Figure 6: 

Comparison of differen t fitting codes, SPS models, and priors on derived stellar masses (from 
iMoustakas et al.ll2012f) . The fiducial m asses are based on fitti ng SDSS and GALEX photometry of 
z ~ galaxies using the iSEDf it code llMoustakas et al.ll20T2ft . with SSPs from FSPS 
jConrov. Gunn fc WhitdliuCih . including dust attenuation, a range in metallicities, and SFHs 
with both smooth and bursty components. The left panel compares stellar mass catalogs produced 
by different groups/codes. K-correct and MPA/JHU-DR7 are based on SDSS photometry; 
MPA/JHU-DR4 is based on SDSS spectral indices, and Salim+07 is based on SDSS and GALEX 
photometry. The middle panel shows the effect of different SPS models (i.e., different SSPs), and 
the right panel shows the effect of varying the priors on the model library. The mean systematic 
differences between mass estimators is less than ±0.2 dex. Figure courtesy of J. Moustakas. 



see Section 13.21 Stellar masses appear to be the most robust parameter es t imated from 
SEP fitting (e.g 



20091 ; iMuzzin et al 



Papovic 
" ' 20091 : 



l. Dickinson fc Ferguson! l200ll ; Ishaplev et"al1l200ll ; IWuvts et al 
Lee et al.l l2009bl ). The reason for this appears to not be fully 



understood, although it probably is at least partly due to the fact that the dust reddening 
vector is approximate ly parallel to SFH and metallicity variations in color-M/L diagrams 
( Bell fc de Jonel 20011 ) . The choice of the SFH, in particular whether it is rising, declining, 
or bursty, can significantly change the best-fit s tellar mass, by perhaps as much as 0.6 dex in 
extreme cases |Pforr. Maraston fc Toninj|2012r >. A general rule of thumb is that M/L ratios 
esti mated via simple SFHs (or single-age models) will be lower l i mits to the true M/L ra- 



tios ( Papovich. Dickinson fc Ferguson 2001 ; Shaplev et al.l 2001 ; Trager. Faber fc Dressier 



200S 



ii lGraves fc Faberll2010l : IPforr. Maraston fc Toninil I2OI2I ). This is a consequence of the 



fact that young stars outshine older ones, making it relatively easy to 'hide' old stellar 
populations in galaxies with a large number of young stars. This also explains why color- 
M/L relations are so uncertain for very blue colors. This general rule is not universal: the 
modeling of certain g alaxy types may lead to systematic biases in the opposite direction 
( Gallazzi fc Bellll2009l ). Partly for these reasons, stellar masses esti mated for quiescent sys- 
tems should be more reliable than for star-forming ones. Relatedly, IZibetti. Chariot fc Rix 
( 20091 ) demonstrated that stellar masses estimated for individual nearby galaxies based on 
a pixel-by-pixel analysis of colors are generally larger than those estimated from integrated 
colors. The differences are typically ~ 0.05 — 0.15 dex, depending on galaxy type. This 
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should not be surprising, especially for systems with star-forming disks and quiescent bulges, 
since the M/L ratio estimated from the integrated colors w ill tend to be biased low by the 
more luminous young stars in the star-forming component. IWuvts et al.l (J2012J) performed 
a similar analysis on galaxies at 0.5 < z < 2.5 with WFC3 Hubble Space Telescope (HST) 
photometry and found no systematic difference in stellar masses derived from resolved vs. 
integrated light, although the fact that their galaxies were at high redshift means that they 
were probing spatial scales an order of magnitude larger than in Zibetti et al. 

There has been some confusion in the literature regarding the importance of restframe 
NIR photometry for estimating stellar masses. First, as indicated in Figure [5] for smoothly 
varying SFHs, photometry at and redward of the V— band is sensitive to the same light- 
weig hted age, and so red der bands do not provide stronger constraints on the mean stellar 



age. I Taylor et al.l (|201ll ) analyzed mock galaxies and concluded that the addition of NIR 
data did not yield more accurate masses. Taylor et al. also found that different SPS models 
produced good agreement in derived properties when NIR data was excluded from their 
fits, but poor agreement when NIR was included. These authors also found much larger 
residuals in their SED fits when NIR data were included, suggesting that the models are 
still poorly calibrated in this regime. As discussed in Section 15.21 the NIR is at present 
probably most useful for constraining metallicities (within the context of a particular SPS 
model), and so NIR data may be useful in cases where there is a degeneracy between M/L 
and Z. In general however stellar mass estimates do not appear to be strongly improved 
with the addition of NIR data, at least with currently available models. Exceptions to this 
rule may be made for gal axies with very high d ust o pacities. 



As first emphasized bv lBell fc de Jond (|200lh and iBell et aD |2003l 'l. some of the largest 



uncertainties in derived M/L ratios stem from uncertainties in the assumed SFHs, in par- 
ticular the presence of bursty SF episodes. The consideration of the Balmer lines with other 
age and metallicity-sensitive features, such as the 4000 A break (D n 4000), can constrain the 
burstiness of the SFH. Optical spectra therefore offers the possibility of providing stronger 
constraints on the M/L ratio. Kauffmann et al. (2003) modeled the H<5 and D n 4000 spec- 
tral features measured from SDSS spectra in order to constrain SFHs and M/L ratios. They 
obtained 95% confidence limits on stellar masses of ~ 0.2 dex and ~ 0.3 dex for quiescent 
and star-forming galaxies, respectively. Again, these are statistical uncertainties because 
the underlying SPS m odel and oth e r asp ects such as the adopted parameterization of the 



SFH were held fixed. IChen et al] (|2012f ) employed principle component analysis (PCA) 



to model the optical spectra of massive galaxies in SDSS. These authors found systematic 
uncertainties in the recovered stellar masses of order ~ 0.1 dex depending on the assumed 
metallicity, dust model, and SFH. They also demonstrated that their PCA technique was 
capable of measuring parameters at much lower S/N than direct fitting to selected spectral 
indices. They determined statistical uncertainties on their mass measurements to be ~ 0.2 
dex, which they showed to be comparable to the formal errors estimated from modeling 
broadband photometry. This result is in agreement with Gallazzi fc Belli ( 20091 ) . who ar- 
gued that for galaxies with simple SFHs and lacking dust, the uncertainty on the derived 
masses are not much larger when using color-based estimators compared to spectroscopic- 
based estimators. With regards to estimating M/L, the real value of spectra appears to be 
restricted to galaxies with unusual SFHs. 

The comparison between spectroscopically-based and photometrically-based stellar masses 
is informative because the approaches suffer from different, though by no means orthog- 
onal systematics. The most obvious difference is with regards to dust — s pectroscopic 
masses are much less sensitive to dust attenuation than photometric masses. iDrorv et al 
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(hood) compared t heir o wn photometrically-derived masses to the spectroscopic masses from 



Kauffmann et al.l (|2003h . They found an rms scatter of ~ 0.2 dex bet ween the two esti 



mates , with a modest systematic trend that correlated with Ha EW. iBlanton fc Roweis 
( 20071 ) estimated photometric stellar masses via a technique that is similar to PC A except 
that the templates are constrained to be non-negative and are based on an SPS model. 
They compared their derived stellar masses to those of Kauffmann et al. and found good 
agreement, with systematic trends between the two restricted to < 0.2 dex. 

A comparison between stellar masses for the same galaxies esti mated with different SP S 
models, priors, and fitting techniques is shown in Figure [6] (from Moustakas et al. 20121 ). 
The galaxies included in this figure are predominantly 'normal' star-forming and quiescent 
z ~ galaxies. For these galaxies the mean absolute differences between various mass 
estimators is less than 0.2 dex, in agreement with other work on the systematic uncertainties 
in stellar masses of normal galaxies. 



3.2 Uncertainties in M/L due to Stellar Evolution Uncertainties 

Maraston et al. ( 20061 ) was the first to draw attention to the sensitivity of deri ved stellar 



masse s to uncertain stellar evolutionary phases, in particular the TP-AGB. The iMaraston 
(2005) SPS mod el predicts much more luminosity arising from TP-AGB sta rs than in previ- 
ous models (e.g., Fioc fc Rocca-Volmerange] 1997 . Bruzual fc Chariot 2003). At ages where 
the TP-AGB phase is most prominent (~ 3 x 10 8 — 2 x 10 9 yr), Maraston' s model predicts 



M araston et al 



rough ly a factor of two more flux at > 1 fim compared to earlier work 
( 20061 ) found that her SPS model im plied a factor of two small er stellar mass, on average, 
compared to masses derived with the Bruzual fc Chariot ( 20031 ) SPS model for galaxies at 
z ~ 2. When the fits excluded the possibility of dust reddening, Maraston's model pro- 
vided a better fit to the optical-NIR SEDs than previous models, whereas when dust was 
allowed, the quality of the fits became indistinguishable, although the offsets in best-fit 
masses remained. 

Subsequent work h as largely confirmed the sensitivity of estimated stellar masses to the 
adopted SPS model JWuvts et al.ll2007l. Kannappan fc Gawisenl2007l ICimatti et al] 120081 . 



Muzzin et al.ll20Qgl . lLonghetti fc Saraccd 



2009, 



Conrov. Gunn 



White! 120091 ). Essentially 
all work on this topic has focused on comparing Maraston's model to Bruzual & Chariot's, 
where the difference in TP-AGB treatment it probably the most significant, though not the 
only difference (other differences include different RGB temperatures and different treat- 
ments for core convective overshooting). Many authors have found a maximum factor of 
~ 2 — 3 difference in derived M/L ratios be t ween the two models, especially when NIR 



data were included. iKannappan fc Gawiserl (|2007l ) showed that the differences between 



Maraston's an d Bruzual fc Chariot's models are relatively modest (< 1.3) when NIR data 
are excluded. Conrov. Gunn fc White ( 20091 ) constructed a new SPS model in which the 
luminosity contribution from the TP-AGB phase could be arbitrarily varied. These authors 
isolated the importance of the TP-AGB phase and confirmed previous work indicating that 
the adopted weight given to this phase in the models can have a large modulating effect on 
the stellar mass. 

The contribution of TP-AGB stars to the integrated light peaks at ~ 3 x 10 s — 2 x 10 9 
yr, depending on metallicity, and so the importance of thi s phase to mode l ing S EDs will 
depend on the SFH of the galaxy. It was for this reason that IMaraston et al. I (|2006l ) focused 
their efforts on quiescent galaxies at z ~ 2; at this epoch even a quiescent galaxy will have a 
typical stellar age not older than several Gyr. In addition to high-redshift quiescent galaxies, 
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Lancon et al. ( 19991 ) suggested that post-starburst galaxies should harbor large numbers of 



TP-AGB stars because their optical spectra show strong Balmer lines, indicative o f larg e 
numbers of 10 s — 10 9 yr old stars. This f act was exploited by IConrov fc Gunnl (201CJ), 



Kriek et al] (fcoid l. and lZibetti et al] (|2012h to constrain the TP-AGB contribution to the 



integrated light. These authors all find evidence for a low contribution from TP-AGB stars, 
and in particular they argue that Maraston's models predict significantly too much flux in 
the NIR for these objects. 

The situation with TP-AGB stars is extremely complex owing to the fact that this phase is 
so sensitive to age and metallicity. Moreover, using galaxy SEDs to constrain the importance 
of this phase is difficult because many parameters must be simultaneously constrained 
(meta llicity, SFH, dust, etc .). A case in point is the analysis of high-rcdshift galaxies, 
where Maraston et al. ( 20061 ) found large differences in the quality of the fit of Bruzual fc 
Chariot models depending on whether or not dust was included in the fits (as an aside, if 
these galaxies do have copious numbers of TP-AGB stars then one may expect them to also 
contain dust, since these stars are believed to be efficient dust factories). It is also worth 
stressing that the conclusion to the TP-AGB controversy may not be an 'either-or' situation 
in the sense that some models may perform better for some ranges in age and metallicity 
while other models may perform better in different regions of parameter space. Ultimately 
of course we desire models that perform equally well over the full range of parameter space, 
and this requires a continual evolution and improvement of SPS models. 

Of course, there are other aspects of stellar evolution that are poorly constrained, includ- 
ing convective overshooting, blue stragglers and the HB, and even the temperature of the 
RGB. The propagation of these uncertainties into derive d properties such a s M/ L ratios 
is only just beginning (e.g JConrov. Gunn fc Whitelbood ). iMelbourne et al.1 (|2012h has for 
example demonstrated that the latest Padova isochrones fail to capture the observed flux 
originating from massive core He burning stars as observed in nearby galaxies. These stars 
are very luminous and can dominate the NIR flux of young stellar populations (< 300 Myr). 

Finally, it is worth stressing that the effect of TP-AGB stars is limited to galaxies of 
a particular type, in particular those that are dominated by stars with ages in the range 
~ 3 x 10 8 — 2 x 10 9 yr. Examples include post-starburst galaxies, which are rare at most 
epochs, or quiescent galaxies at z > 2. For typical galaxies at z ~ 0, the treatment of this 
phase seems to be of little relevance for estimating stellar masses, at least on average, as is 
evident in Figure [6] 



3.3 Stellar Masses at High Redshift 



As first emphasized bv iPapovich. Dickinson fc Ferguson! ([20011 ). galaxies with high SFRs 
can in principle contain a large population of 'hidden' older stars, since these stars have 
high M/L ratios. Papovich et al. found that the data allowed for significantly larger stellar 
masses when two-component SFH models (young+old) were compared to their fiducial 
single component models (a median difference of a factor of ~ 3 with extreme cases differing 
by an order of magnitude) . This issue becomes more severe at higher redshifts because the 
typical SFRs are higher and because of the increasing possibility of rising SFHs (to be 
discussed in detail in Section [4}. This makes the analysis of high redshift galaxy SEDs 
much more complicated. 

However, at the highest redshifts (z > 6), the analysis of SEDs may actually become 
simpler. At z = 8 the age of the universe is only « 640 Myr, and if one presumes that 
galaxy formation commences after z m 20, then even the oldest stars at z — 8 will be no 
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more than « 460 Myr old. The oldest possible main sequence turnoff stars will thus be A 
type stars. Put another way, the M/L ratio in the blue (U through V bands) at 13 Gyr 
is 4.5 — 15 times higher for an instantaneous burst of SF compared to a constant SFH, 
whereas the M/L ratio for these two SFHs differs by the more modest 2.5 — 5 at 500 Myr. 
So there is some expectation that 'hidden mass' will be less hidden when modeling SEDs 
at the highest redshifts . This expectation appears to be borne out by detailed modeling. 



Finkelstein et al. (2010) analyzed the SEDs of z ~ 7 — 8 galaxies and concluded that even 



allowing for the possibility of extreme amounts of old stars (90% by mass) , the best-fit stellar 
masses increa sed by no more than a f actor of two compared to a fiducial single-component 



SFH (see also ICurtis-Lake et al.ll2012r ) 



At even higher redshifts [z > 10), the analysis of galaxy SEDs may become even simpler, 
at least with regards to estimating stellar masses. At sufficiently high redshifts the age 
of the oldest possible stars will eventually become comparable to the SF timescale probed 
by the restframe UV (~ 10 s yr at ~ 2500 A; see Figure [5]). The measured SFR will thus 
be averaged over roughly the entire age of the galaxy, implying that a reasonable estimate 
of the stellar mass can be obtained by multiplying the UV-derived SFR with the age of 
the universe. This may in fact explain the unex pectedly strong c orrela tion between stellar 



mass and UV luminosity at z > 4 reported by iGonzalez et al.l (|201lf ). Of course, at the 
highest redshifts other difficulties arise, including accounting for the effects of very high 
EW emission lines and the reliability of the underlying models at very low metallicities. 



3.4 Summary 

M/L ratios for most galaxies with normal SEDs are probably accurate at the ~ 0.3 dex level, 
for a fixed IMF, with the majority of the uncertainty dominated by systematics (depending 
on the type of data used). Galaxies with light-weighted ages in the range of ~ 0.1 — 1 
Gyr will have more uncertain M/L ratios, with errors as high as factors of several, owing 
to the uncertain effect of TP-AGB stars and perhaps other stars in advanced evolutionary 
phases. Galaxies that have very young light-weighted ages (e.g., high-redshift and starburst 
galaxies) will also have very uncertain M/L ratios because of the difficulty in constraining 
their past SFHs. Uncertainties due to dust seem to be subdominant to other uncertainties, 
at least for standard reddening laws and modest amounts of total attenuation. Essentially 
all stellar masses are subject to an overall normalization offset due both to the IMF and the 
uncertain contribution from stellar remnants. Owing to the lingering model uncertainties 
in the NIR (e.g., from TP-AGB and cool core He burning stars), it may be advisable to 
de rive physical param eters by modeling the restframe UV-optical (see e.g., the discussion 
in iTavlor et allboilh . 



4 STAR FORMATION RATES, HISTORIES, & STELLAR AGES 
4.1 Measuring SFRs 

The practice of measuring SFRs from mono chromatic indicators (and combinations ther eof) 
is well established (see e.g., the reviews of lKennicuttlll998l , iKennicutt fc Evansll2012h . A 



wide variety of relations have been established between SFRs and the UV, Ha, Pa, 24 
luminosity, total IR luminosity, radio continuum luminosity, and even the X-ray flux. As 
they are not the focus of the present review, only a few brief remarks will be given regarding 
monochromatic SFR indicators. Each indicator is sensitive to a different SFR timescale, 
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Figure 7: 

Comparison between SFRs estimated from e mission lines ( "Ha " ; iBrinchmann et aljliool) and 
from modeling the UV-optical SEDs ( "UV" ; ISalim et alj|200^ of star- forming galaxies at z ~ 0. 
The right panel shows the difference between SFR estimators as a function of stellar mass. There 
is a systematic trend between the two estimators, which may be due to modeling differences 
and/or the uncert ainty in extrapolating the fiber-based Ha SFRs to the total galaxy. From 
ISalim et all l|2007l l. 



with ionizing radiation (as probed by e.g., the recombination lines) being sensitive to the 
shortest timescale, ~ 10 6 yr, and the UV and total IR being sensitive to longer timescales, 
e.g., ~ 10 8 yr. In detail, the UV and IR are sensitive to a range of SF timescales depending 
on wavelength and SFH. In the IR this sensitivity arises because shorter wavelengths are 
probing dust heated mostly by H II regions, while FIR photons are sensitive to dust heated 
both by old and young stars. Amongst the recombination lines, 'redder is better' since 
dust attenuation will be less severe at longer wavelengths, and so Pa should be preferred 
over Ha. For reference, an attenuation of 1 mag in the V— band implies an attenuation of 
only 0.15 mag at Pa, for standard attenuation curves. The optical and UV SFR indicators 
require correction for dust when used on their own. For this reason, mixed SFR indicators 
have become popular, usually via the comb ination of UV and IR or Ha and IR indicators 



( Calzetti et al.ll2007l . iKennicutt et alj|2007f ) 



Estimating SFRs from SED fitting is very challenging for several reasons: (1) the age- 
dust-metallicity degeneracy makes it difficult to reliably measure ages and hence SFRs 
unless high quality data are available (see Section (2) the choice of model priors on 
the dust model and SFH library imposes often severe biases on the resulting ages and 
SFRs. These difficulties were recogniz ed early on in the study of high-redshift galaxies 
(Papovich, Dickinson & Ferguson 2001) and are still relevant today. 



Brinchmann et al.l ( 20041 ) measured SFRs for SDSS galaxies by modeling a suite of optical 
emission lines. Their model included variation in the total metallicity, dust attenuation, 
ionization parameter, and dust-to-metals ratio. As the SDSS is a fiber-based survey, these 
authors had to make a correction for the fact that the spectra sample only a fraction of the 
galaxy light. In addition, they modeled the st ellar continu u m in order to subtract off any 



absorption features under the emission lines. ISalim et al.l (|2007f ) analyzed a comparably 



sized sample of z ~ galaxies, focusing on UV photometry from the Galaxy Evolution 
Explorer (GALEX) and optical photometry from SDSS. These authors fit the data to SPS 
models with a range of SFHs (including exponentially declining SFHs and superimposed 
starbursts), metallicities, and dust attenuation. They presented a direct comparison be- 
tween their SED-based SFRs and the emission line-based SFRs from Brinchmann et al. 
The result is shown in Figure [7] While the overall agreement is good there is a systematic 
trend with stellar mass resulting in offsets of ~ 0.2 dex at the low and high mass ends. 
Salim et al. explored the origin of this offset and concluded that the different derived dust 
attenuation optical depths between the two methods were the source of the discrepancy. 
They further speculated that in most cases it was the emission line-based dust attenuation 
values that were more accurate, and therefore also the SFRs. Salim et al. also demon- 
strated that for galaxies with low SFRs and/or galaxies with emission line ratios indicative 
of AGN activity, the UV-based SFRs tended to be more reliable than those estimated by 
Brinchmann et al. These results suggest that SFRs based on modeling UV-optical SEDs 
c arry systemat i c unc ertainties at the < 0.3 dex level. 



Wuvts et alj (|201ll ) compared the SFRs of galaxies at z ~ 3 computed with three different 



techniques: Ha-based SFRs, SED-fitting based SFRs, and a mixed SFR indicator based on 
the UV and IR emission. For SFRs < lOOMoyr -1 , these authors found that the three 
indicators agreed reasonably well (no systematic offsets were detected), provided that extra 
attenuation was included toward H II regions when modeling the Ha-based SFRs. They also 
found that at the highest SFRs > 100 Mq yr" 1 the S ED-based values ten ded to underpredict 
the SFR, compared to the mixed indicator (see also Santini et al. 20091 ) . They argued that 



this was due to their assumption of a uniform dust screen, and at high column density this 
becomes an increasingly poor representation of reality. 

At high redshift the effect of nebular emission lines becomes an additional source of 
systematic uncertainty in the modeling, and can significantly affect the best-fit SFRs. The 
impact of emission lines is larger at higher redshift because of three effects: (1) specific 
SFRs are higher; (2) metallicities are lower; and (3) the redshifting of the spectrum implies 
that a line of given EW comprises a larger fraction of the flux through a filter. In Figure 
[8] the specific SFR for galaxies with a mass of 10 9 ' 5 M Q is shown as a function of redshift 
for various choices regarding the contribution of emission lines to the observed broadband 
fluxes. Evidently the inclusion or not of emission lines in the modeling can change the 
derived specific SFRs in a r edshift-dependent manner, by a facto r of roughly two in the 
most extreme cases (see also iPacifici et al.ll2012l . lLabbe et al]|2012h . 



4.2 Light- Weighted Mean Ages 

The mean stellar ages of galaxies are another important probe of the SFH. Of course, 
mass-weighted ages offer a direct probe of the integrated SFH, while light-weight ed ages 
are more directly measurable from SEDs. Mass-weighted ages will always be larger than 
light-weighted ones, owing to the fact that young stars outshine older stars. As discussed 
further in the next section, the mass-weighted ages of actively star-forming galaxies turn 
out to be highly sensitive to the assumed model SFH used in the SED fitting. Due to the 
fact that light-weighting is highly non-linear, light-weighted ages computed from spatially 
resolved photometr y will tend to be old er than ages computed from integrated photometry, 
as demonstrated bv lWuvts et al.l (|201ll ). 

The mean ages of quiescent galaxies have been estimated largely via the modeling of 
spectral features, including D n 4000 and the Lick indices. A widely adopted procedure is 
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Figure 8: 

Specific SFRs for galaxies with stellar masses of 10 9 ' 5 Mq as function of redshift. The black and 
green symbols reflect different assumptions regarding the SFH ('REF', 'RIS', 'DEC for constant, 
rising, and declining) a nd the correction for emission lines ('NEB'). Error bars reflect 68% 
confidence limits. From lde Barros. Schaerer &: Starlj l|2012f ) 



to fit single-age SSPs to the Lick indices in order to estimate ages. In this approach it is 
the hydrogen Balmer lines that provide the primary constraint on the age, as these lines 
are strongest in A type stars and thus for populations with ages of ~ 10 s — 10 9 yr. The 
derived ages are sometimes referred to as 'SSP-equivalent ages' in order to highlight the 
fact that no attempt is made to model composite stellar populations. A number of authors 
have demonstrated t hat these derived ages are lower limits to t he true mass-weighted and 



light- weighted ages (ITrager et all 2000a! lierra fc Trageil |2007|. Frager fc Somervilk |2009| . 



MacArthur. Gonzalez fc Courteaul 12009V As emphasized by ITrager fc Somervillel (|2009r i 



this fact implies that any observed trend of SSP-equivalent ages with other parameters such 
as velocity dispersion will be stronger than the true underlying mass-weighted trend. The 
SSP-equivalent ages are younger even than ligh t-weighted ages largely beca use of the sensi- 
tivity of the Balmer lines to ~ 1 Gyr old stars. ITrager fc Somervillel i|2009ri concluded that 
SSP-equivalent ages are mostly measuring residual star formation over ~ 1 Gyr timescales 



Ferreras et al 



The se conclusions are not restricted to the analysis of spectral indices 
( 20091 ) estimated mean ages of early-type galaxies at z ~ 1 by fitting models to HST grism 
spectra. The ages estimated from single-age models showed almost no correlation with 



2!) 



ages estimated with composite models. Likewise, MacArthur. Gonzalez fc Courteau 1 20091 ) 
performed full spectrum fitting on a sample of nearby galaxies and found only a modest 
correlation between ages derived from single-age models and light-weighted ages derived 
from composite models. 

It is widely believed that the H/3 index is the id eal age indicator because it is largely 
insensitive to metallicity and ab undance patterns (Worthev 1994. Frripicco fc Bell 1995 . 
KOTn i ^arastori_&_TtornasJ_ 2005 \. H owever, as noted by WorthevT i 1994 ) and later by 



de Freitas Pacheco fc Barbuv (1995), Maraston fc Thomas 1 2000h . and Lee. Yoon fe Leel 



(2000), the Balmer lines are also sensitive to the presence of blue HB stars, and there will 
therefore be some ambiguity in determi ning ages from H/3 without additional constraints 
on the presence of such stars. However. iTrager et alj (|2005h demonstrated that, while the 
inferred ages of the oldest stellar populations (> 10 Gyr) are susceptible to uncertainties of 
2 — 5 Gyr due to the presence of blue HB stars, the ages of intermediate-age populations 
(~ 1 — 10 Gyr) are rather less affected by blue HB stars. The reason for this is that 
the Balmer lines increase so strongly with decreasing age that very large quantities of 
blue HB stars would be required to be pre sent in old ste ll ar po pulations to reproduce 
the high observed H/3 EWs. As an example, iThomas et al. |2005h found that they could 
reproduce the highest H/3 EWs in their sample with a model in which 50% of the old 
metal-rich stars had blue HB stars. In this two-component model, the component with 
blue HB stars needed to have a metallicity 0.2 — 0.4 dex higher than the other component. 



Percival fc Salaria (|2011f ) also found that models with very extreme HB populations and 
old ages could masquerade as a relatively young population with ages as young as ~ 3 
Gyr, based on H/3 EWs. There is no compelling evidence supporting the existence of a 
population of blue HB stars in metal-rich quiescent galaxies large enough to explain the full 
range of H/3 EWs. However, it is reasonable to assume that there will be a population of 
blue HB stars in such systems, and they will affect the derived stellar ages. 

The higher order Balmer lines are promising because they will be even more sensitive 
to warm/hot stars than H/3. In principle then, the consideration of multiple Balmer 
lines should provide constraints on multiple age components, and /or al l ow for the sep- 
arati on of age and blue HB effects (|Schiavon et al.l |2004| . ISchiavonl 120071 . ISerra fc Trager 
2007). However, the potential power of the higher-order Balmer lines is tempered by the 



fact that these features are more sensitive to metallicity and abun d ance ratios than H/3 
( Thomas. Maraston fc Kornlbool iKorn. Maraston fc Thomas! I200H ). ISchiavonl |2007h an- 
alyzed the H/3, H7, and H5 indices measured from stacked SDSS early-type galaxy spectra 
with a variety of composite models including varying blue HB stars, blue stragglers, and 
multi-age populations. He concluded that small fractions of young/intermediate-age stellar 
populations, perhaps in the form of extended SFHs, was the most likely scenario to ex- 
plain the observed trend s. This propo s al is si milar, though not identical to the 'frosting' of 
young stars proposed by Trager et al. ( 2000al ) in order to explain the observed correlations 
between velocity dispersion, age, metallicity, and abundance pattern. 

There exist ot her proposed methods for constraining multi-age components based on 
spect ral indices 
line 



Rose 



Leonardi fc Rose 



19841) proposed a Ca 11 index that is sensitive to the He Balmer 
1996T ) demonstrated that this index in conjunction with H5 can 
break the 'age-str ength degeneracy', referr i ng to the age and strength of a recent burst 
of star formation. Smith. Lucev fc Hudson ( 20091 ) used this diagnostic to argue against a 
frosting of small fractions of young stars (< 1 G yr) to explain the young apparent ages for 
quiescent galaxies in the Shapley supercluster. iKauffmann et al. (2003) considered both 
the H<5 and D n 4000 spectral features in an attempt to measure two-component SFHs (in 
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particular, a burst of SF and an underlying r model SFH). Galaxies with a sizable mass 
fraction formed in a recent burst of SF are significant outliers in the }iS- D„4000 plane, and 
so str ong constraints on the burst fraction can be obtained in this way. Percival fc Salarisl 
( 201 In argued that the Ca n Rose index and the Mg II feature near 2800 A could also be 
u sed to separate age fr om b lue HB effects. 

Renzini fc Buzzoni ( 19861 ) suggested that the onset of the AGB in ~ 0.1 — 1 Gyr old 



po pulations could b e used to age-date galaxies in the NIR. Subsequent modeling efforts 
by iMarastonl (|2005l ) confirmed that NIR spectral features attributable to oxygen-rich and 
carbon-rich AGB stars vary strongly over the age range ~ 0.1 — 1 Gyr, and may be used as 
an independent constraint on the presence of intermediate-age populations in galaxies. In 
principle, the very strong sensitivity of AGB evolution to mass (age) and metallicity suggests 
that they should be good clocks. However, their sensitivity to mass and metallicity is due 
to poorly understood physics, including co nvection and mass-loss, and so the calibration 



of an AGB-based clock is quite uncertain. iMiner. Rose fc Cecill (|201ll ) attempted to use 



NIR spectral features to place constraints on the ages of M32 and NGC 5102 based on 
Maraston's models. They found broad qualitative agreement between NIR-based ages and 
CMD-based ages. Future quantitative modeling of the optical-NIR spectra will be required 
to assess the utility of AGB-sensitive spectral indices as age indicators. 

Finally, it is also worth emphasizing that it is fundamentally very difficult to measure ages 
for systems older than > 9 Gyr because evolution in the isochrones is so slow at late times. 
For example, at solar metallicity between 9 — 13.5 Gyr, the RGB changes by only « 50 K 
while the main sequence turnoff point changes by « 230 K. Very accurate models and very 
high quality data are required to be able to distinguish such small temperature variations. 
In theory the UV should be much more sensitive to age at late times due to the onset 
of an extended (blue) HB (|Yi et alJll999h . In the context of a particular model, the UV 
colors become very sensitive to age for ages > 7 Gyr. The problem with this chronometer is 
that it is very sensitive to the uncertain physical inputs, such as the mass-loss prescription 
and the helium abundance, an d also to the underlying distribution of stellar metallicities 
( Yi. Demaraue fc Oemlerl ll9970. These uncertainties will render absolute ages unreliable, 
although relative ages (or even simply the rank-ordering of galaxies by age) are probably 
more robust. 



4.3 Which Model SFHs to Use? 



Papovich. Dickinson fc Ferguson! (|200ll ) recognized that SED modeling suffered from a fun- 
damental limitation, namely that the young stars, being so bright, tend to outshine the 
older, less luminous stars. This outshining effect was known t o result in underestimated stel- 
lar masses when single-component SFH models are used (e.g., Papovich. Dickinson fc Ferguson! 



2001 



Daddi et al 



2004; Ishaplev et all I2OO5I ; iLee et all l2009bl ; IPforr. Maraston fc Tonini 



20121 ) . b ut its effect on rec overed SFRs and ages remained relatively less well-explored until 
recently. Lee et al. ( 2009bl ) analyzed the SEDs of high-redshift mock galaxies and concluded 
that the use of single-component exponentially-decreasing SFHs (r models) resulted in sub- 
stantially underestimated SFRs and overestimated mean ages (by factors of two for both 
parameters) . They attributed this shortcoming to the fact that the mock galaxies tended to 
have r i sing SFHs in contrast t o the declining mode l SF Hs. Subsequent work by 



Lee et al 



( 20ld l lMa~raston et all <|2010| ). IWuvts et all I (|201lh . and IPforr, Maraston fc Toninil <|201^ T 
largely confirmed these findings, and further concluded that models with rising SFHs tended 
to provide a better fit to observed high redshift SEDs and produced SFRs in better agree- 
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ment with other indicators. Relatedly, the choice of the allowed range of each SFH pa- 
rameter (i.e., the prior) can signi ficantly affect the inferred best-fit SFRs and ages (e.g., 



Lee et al 



20ld . lWuvts et alj|201ll ). 



A specific example of the effect of the adopted SFH parameterization on SFRs is shown 
in Figure [8] In this figure the specific SFR is shown as a function of redshift at a reference 
stellar mass of 10 9 ' 5 M Q . At high redshift the scatter between the "RIS", "DEC", and 
"REF" models is entirely due to the choice of the assumed SFH (rising, declining, or 
constant, respectively). At the highest redshifts these different choices result in a change in 
specific SFRs of more than an order of magnitude. 

There has been some diff erence in opinio n rega rding which functional form one should 
adopt for the rising SFHs. iMaraston et all (hold). andlPforr. Maraston fe Toninil (|2012l ) 
advocated exponentially increasing SFHs while lLee et ail (|201ol ) advocated delayed r models 
of the form SFRoc t e _t,/r . The data do not presently favor one functional form over the 
other. The basic conclusion to draw from these analyses is that the model SFH library 
must be sufficiently diverse to allow for a wide range in SFH types. Model SFHs that are 
too restricted will inevitably lead to biased results. Of course, the criterion for 'sufficiently 
diverse' SFHs will depend on the galaxy type in question and has not been thoroughly 
examined in a general sense. At present, all that can be stated with confidence is that for 
certain SED types, especially highly star-forming galaxies, great care must be taken when 
constructing the model library of SFHs. 

Another option is to use a library of SFHs drawn from s emi-analytic or hydrodynamic 
model s of galaxy formation. This technique was introduced bv lFinlator. Dave fc Qppenheimer 
(2007|), who analyzed SEDs of six high redshift galaxies with SFHs drawn from their hydro- 
dynamic simulations. The advantage of this approach is that the model SFH library is more 
likely to contain realistic SFHs ( with rising and fallin g SFHs, multiple bursts, etc.) than 
simple analytic models. Recently, iPacifici et al.l (|2012l ) explored a similar technique using 
SFHs from a semi-analytic model rather than a hydrodynamic simulation. The advantage 
of using a semi-analytic model is that a much larger range of SFHs can be computed in a 
given amount of computational time. As with any other SFH library, the results can be 
sensitive to the choice of priors. In this case, the 'prior' is in some sense the reliability of 
the model in producing a realistic range of SFHs. A possible disadvantage of this approach 
is that it can be much more difficult to quantify the influence of this prior on the derived 
physical properties than when using analytic SFH libraries. 



4.4 Non-Parametric Star Formation Histories 



When modeling SEDs one typically must adopt some parameterization of the SFH, and 
it is this step that clearly introduces a number of poorly quantified systematics. An ap- 
pealing solution to this problem is to model SEDs with non-pa rametric SFHs. This is pre- 



cisely the approach ta ken by several groups, inclu ding MOPED ( Heavens. Jimenez fc Lahav 



2000 ), STARLIGHT dCid Fernandes et alj|2005l) STEC MAP (lOcvirk et al.ll2006D. VESPA 
( Toieiro et al.ll2007l ). ULvSS jKoleva et al.l l2009h . and iMacArthur. Gonzalez fc Courteau 
( 2009)) . In the simplest implementation one specifies a fixed set of age bins and then fits for 
the fraction of mass formed within each bin. More sophisticat ed algorithms utiliz e adap- 



tive age binning that depends on the SED type and S/N (e.g.. iToieiro et al.ll2007r ). In all 
cases one is engaging in full spectrum fitting (in the restframe optical; for an example of 
the technique see Figure [9]). In principle one could use such codes to fit broadband SEDs, 
but the fits would likely be seriously under-constrained. These methods have been tested 
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Figure 9: 

Example fits to the optical spectrum of a quiescent and a star-forming galaxy with VESPA 
(Toic rro et a,l.|[2009t) ■ The SFHs are estimated non-parametrically with an age binning that is 
adaptive depending on the S/N and information content of the spectrum. The upper panels show 
the best-fits (red line is the model, black line is the data; emission lines have been excluded from 
the fits); lower panels show the recovered SFHs plotted as the mass fraction formed within each 
age bin. Figure courtesy of R. Tojciro. 



in a variety of ways, and it appears that they are capable of recovering complex SFHs 
remarkably well, at least when the mock galaxies are built with the same SPS mod els as 
used in the fitting routines (jOcvirk et al.ll200d . lToieiro et al.ll2007l . lKoleva et alj|200gh . Not 
surprisingly, the requirements on the data quality are demanding: a wide wavelength cov- 
erage, high spectral resolution, and high S/N (typ i cally > 50/A) are required in order to 



robustly recover complex SFHs ( Ocvirk et al. 2006I . Toieiro et al. 2007 ). In addition, there 



are clear systematics in the recovered SFHs due to different SP S models, particularly for 

~ 20091 



the SFHs at ages of ~ 0.1 — 1 Gyr (I Toieiro et al 



20111 ). This is perhaps related 



to the different treatment of AGB stars and/or core convective overshooting amongst SPS 
models. Moreover, current implementations of non-parametric SFH recovery utilize SPS 
models that only allow variation in metallicity and age — the abundance pattern is fixed to 
the solar value. When fitting to spectra this practice can introduce additional systematics 
since many of the features that are being fit are sensitive to the detailed abundance pattern 
of the system. 

There have been few comparisons between the SFHs derived via non-parametric methods 
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Figure 10: 



Comparison between SFHs estimated from integrated light SE Ds to SFHs est i mated from CMDs. 
The SED-based SFHs are based on SDSS spectr a, as derived inlToieiro et al ] 1120091 1. and binned 
by stellar mass. The CMD-based SFHs are from lWeisz et at based on data from the 

ANGST survey and are for low-mass dwarf irregulars. It is striking that the CMD-based SFHs fo r 
low-m ass dwarfs look closer to the SED-based SFHs for the most massive galaxies. From llTeitnerl 
(2012h . 



with more conventional techniques. Toieiro et al. d2009h compare d SFRs estimated with 
VESPA to emission line-based SFRs from lBrinchmann et al.l (|2004h . No obvious discrepan- 
cies were found, but further detailed analyses would be desirable. Moreover, there appears 
to be no published work comparing the mean stellar ages (or some other moment of the 
SFH) between these various techniques. 



SFHs from CMDs versus SFHs from Integrated Light Ideally one would like an 
independent constraint on the SFH to compare against the SFH derived from integrated 
light spectra. Unfortunately no direct tests of this sort have been performed to date. The 
most robust method for determing SFHs is via the modeling of resolved color magnitude 
diagrams of the stars within galax ies. This technique is well-developed, and has recently 
been applied bv lWeisz et al.l (|201ll l to a large sample of dwarf gala xies in the local volum e 
compiled by t he ACS Nearby Galaxy Treas ury Survey (ANGS T; [Palcanton et al"ll2009T l. 
Leitner (|2012h compared t he SFHs derived bvlWeisz et alJ ((20 1 In to the SFHs derived from 
integrated light spectra by Toieiro et al. ( 20091 ) . The result is shown in Figure [TO] The dwarf 
galaxies from the ANGST Survey have masses slightly lower than the lowest mass bin in 
the integrated light sample, and it is worth emphasizing that there is no overlap between 
the two samples. And yet, the dwarfs have SFHs comparable to the highest mass bin in 
the integrated light sample. There are many complications in making this comparison, not 
the least of which is the question of whether or not nearby dwarfs are representative of 
the cosmological average. Nonetheless, the difference in the qualitative trend derived from 
CMD-based and integrated light-based SFHs is striking and worthy of further study. 
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4.5 Low-Level Star Formation 



Probing low levels of star formation in galaxies is challenging because at low levels classic 
SF indicators can become contaminated by other processes. At low luminosities, emission 
lines can be influenced by low-level AGN activity and they become difficult to measure on 
top of the strong absorption features in the continuum, dust emission can be due to heating 
from old stars (rather than tracing the formation of new stars), and ultraviolet light can be 
attributed to hot evolved low-mass stars. The former two indicators have been effectively 
abandoned as low-level SF indicators, while the third, UV light, still holds some promise. 

It is beyond the scope of this review to discuss in detail the complications in interpreting 
UV (especially FUV) emission from predominantly quiesce nt systems, but a few remarks will 
serve to round out this section. The reader is referred to lO'Connell |l999h for a thorough 
review, and [3 ( 20081 ) for a recent update. Low-level UV emission from galaxies can be 
due either to young massive stars or to hot evolved low-mass stars. In the latter category 
includes extreme HB stars, post-AGB stars, AGB-Manque, etc. Extreme HB stars are 
known to exist in essentially all low metallicity and some high metallicity globular clusters. 
Post-AGB stars should exist at all metallicities, and AGB-Manque are thought to occur only 
in special circumstances. It is very difficult to predict the specific numbers of these stars as 
a function of metallicity and time, and it is for this reason that low-levels of UV flux are so 
dif ficult to interpret. Th e recent catalog of evolved hot stars in Galactic globular clusters 
bv lSchiavon et alJ (|2012l ) may help to calibrate the models for these advanced evolutionary 
phases. 

■ have been several important developments since the review of lO'ConnelJ l| 19991 ). 
which are mainly due to progress provided by HST and GALEX. It is now clear that galaxies 
on the optically-defined red sequence (e.g., in B — V or g — r c olor space) frequently harbor 
low levels of star formation, at the rate of ~ 0.1 — 1 Mg yr - 1 ( Yi et al. 20051 . Kavirai et al 



2007l . lSalim fe Richl2010l . lSalim et aHl20lllFang et al.ll2012h . although for the most massive 



The implication is that optical SEDs are not 

-l 



galaxies SF levels appear to be even lower. 

sensitive to specific SF at levels of < 10 -11 yr _1 . It is also now quite clear that optically 
defined red sequence galaxies have a wide range of UV colors, from those with low levels of 
star formation to those that are extremely red. It has become common practice to define 
an arbitr ary UV-optical color boundary i n order to separate star forming from quiescent 



galaxies l Yi et al. 2005 . Jeong et al. 20091 ) . Of course, such a practice does not provide an 



answer to the question of how low star formation rates actually are in the reddest galaxies. 
Instead it signals the difficulty in measuring SFRs for such red systems because any UV 
flux becomes easily confused with hot evolved low-mass stars. Further progress can be 
made by directly resolving luminous young stars and star clus ters in nearby galaxies. This 
was precisely the strategy taken by Ford fc Breeman 1 20121 ). who observed four nearby 
elliptical galaxies with HST. These authors argued that the galaxies in their sample have 



current specific SFRs in the range of 10 



10" 



yr 



far below what can be reliably 



probed via SED modeling. UV spectroscopic features may also aide in disti nguishing hot 
evolved low-mass stars from hot upper main sequence stars ( O'Connell 19991 ). 



4.6 Summary 

Measuring the detailed SFH of a galaxy from its observed SED is one of the holy grails of 
SPS modeling. Unfortunately, two unrelated physical effects make this task very difficult 
in practice: (1) high mass stars are so luminous that they dominate the SEDs of star- 
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forming galaxies, easily outshining the more numerous, older, lower mass stars; (2) at 
late times isochrones evolve very little, so that it is difficult in practice to distinguish an 
~ 8 Gyr population from a ~ 12 Gyr one. The first effect implies that light-weighted 
ages will always underestimate the true mass-weighted age. It also implies that the priors 
on the model SFHs can strongly influence the best-fit parameters. Non-parametric SFH 
reconstruction methods offer the best hope of providing unbiased results, but they require 
data of very high quality. As a general rule of thumb, optical SEDs are sensitive to specific 
SFRs as low as 1CT 11 yr _ \ UV fluxes are sensitive to levels as low as 10 12 yr 1 , while 
resolved UV photometry is in principle sensitive to specific SFRs down to at least < 1CP 14 
yr -1 . SED-based SFRs are probably accurate at the factor of two level (at a fixed IMF), 
except for extreme systems where uncertainties can be much larger. There appears to be 
little consensus on the accuracy of light-weighted and mass-weighted stellar ages derived 
from SED fitting. 



5 STELLAR METALLICITIES & ABUNDANCE PATTERNS 



5.1 The Age-Metallicity Degeneracy 

The metallicity, Z, affects the SED in two distinct ways: (1) an increase in metallicity results 
in lower effective temperatures, including a cooler main sequence and giant branch; (2) at 
fixed T c ff, an increase in metallicity results in stronger spectroscopic absorption features 
and generally redder colors. Both of these effects contribute to the overall reddening of 
an SED with increasing metallicity. In addition to the overall metal content, the detailed 
elemental abundance pattern is imprinted on the SED, especially in the strength of spectral 
absorption features. 



Age will also tend to redden the SED through the effect of age on the isochrones. IWorthev 



(|l994f ) studied the degeneracy between age and metallicity in detail using his SPS models 
and concluded that at ages > 5 Gyr neither broadband colors nor most spectral indices were 
able to reliably separate age from metallicity. Worthey introduced his "3/2 rule", whereby 
an increase/decrease in the population age by a factor of three is almost perfectly degenerate 
with an increase/decrease in metallicity by a factor of two. He identified a small number 
of spectral indices that were unusually sensitive to age, including the hydrogen Balmer 
lines, and several indices tha t were unusually sensitive to metallicity, in cluding the Fe4668 
and Fe5270 indices (see also IO'Connei]|l98d iRabinl |1982|, I Rose! 1 19851 , for early attempts 
at separating age and metallicity). Worthey therefore demonstrated conclusively that the 
age-metallicity degeneracy could be broken by jointly considering a handful of carefully 
chosen indices. He also pointed out that the age-metallicity degeneracy is less severe at low 
metallicity, where the turnoff point and giant branch have very different age sensitivity at 
fixed metallicity. 

It has long been appreciated that a clean separation of ag e and metallicity can be achieved 
with Lick indices only for data of very high S/N (> 100/A: Frraeer e t al. 2000b, Kuntschner 
2000). Figure [Til demonstrates that an age-metallicity degeneracy persists when exploiting 



Lick indices at lower S/N. In contrast, full spectral fitting offers a much more robust separa- 
tion of age and metallicity at moderate to low S/N (at least within the context of an underly- 
ing s et of model SSPs; see MacArthur. Gonzalez fc Courteau 20091 ; Sanchez-Blazquez et al. 



2011). This is simply a consequence of the fact that much more information is available in 



the full optical spectrum than can be captured by a limited number of spectral indices. 
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Figure 11: 

Illustration of the age-metallicity degeneracy in spectra of modest S /N. Three techniques are 
compared, the full spectrum fitting technique of lOcvirk et al.l l|2006l . left panel), classic two-index 
fitting (middle panel), and multi-index fitting (right panel). Results are shown for a mock galaxy 
with S/N= 50/A, [M/H]= 0.0, and an age of 10 Gyr (top panels) and 1 Gyr (bottom panels). 
Each point represents the best -fit for a single realization of the mock galaxy spectrum. From 
ISanchc z-Blazquc z et al] l|201ll ). 



5.2 Photometric Metallicities 



The age-metallicity degeneracy discussed in lWorthevI (|l994T ) strictly applies to systems that 
are coeval with ages > 5 Gyr; models at younger ages were not considered by Worthey. In 
subsequent work it became clear that the broadband optical-NIR colors of composite stellar 
populations could separate age and metallicity effects. The physics underlying the idea is 
that blue colors probe the age-sensitive main sequence turnoff point whi le red/NIR colors 



probe the metallicity-sensitive giant branches. Bell fc de Jonj ( 200d ) and MacArthur et al 
(120041 ). amongst others, showed that age and metallicity vectors were nearly orthogonal 
in optical-NIR colors (e.g., B — R vs. R — H). These authors therefore concluded that 
estimates of stellar metallicities were possibl e from photom e tric d ata alone, for systems 
whose light- weighted ages were not too old. Eminian et al. ( 20081 ) showed that galaxies 
grouped by spectroscopically-derived ages and metallicities separate relatively well in g — r 
vs. Y — K color-color space. Notice that this statement does not depend on the predicted 
model colors, thereby providing strong evidence that separation of age and metallicity is at 
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le ast in principle p ossible based on broadband optical-NIR data alone. 

(2003) investigated the model uncertainties associated with separating age and 



Lee et al 



metallicity based on broadband colors. They echoed previous work in that the young 
populations in composite stellar models dominate the light and therefore allow a separation 
of age and metallicity effects when optical-NIR colors are employed. However, they also 
emphasized that different SPS models pr oduce very diffe r ent ag e-metallicity vectors in color- 

The result is that age and 



color space, a point also emphasized bv lEminian et al 



metallicity can be separately constrained within the context of a particular SPS model, but 
the derived ages and metallicities will vary from model to model. The underlying cause 
of this problem is that the NIR light is heavily influenced by AGB stars for populations 
of moderage age, and so different treatments of this uncertain evolutionary phase cause 
substantial differences in the age-metallicity grids. Because of the difficulty in measuring 
metallicities from photometric data, many SPS studies simply treat metallicity as a nuisance 
parameter to be marginalized over. Other studies take the even simpler approach of fixing 
the metallicity to the solar value. The danger in making the latter simplification is that 
assumptions in the assumed metallicity will generally affect other par ameters of interest , 
such as the stellar mass, dust opacity, and star formation rate (see e.g., W uvts et al. 2009, 
Marchesini et al.ll2009l , iMuzzin et allboogl . IPforr. Maraston fc Toninill2012l 'l. 



5.3 Spectroscopic Metallicities 

The more robust approach to estimating stellar metallicities is to employ spectroscopic 
features. Optical spectra of galaxies are rich in atomic and molecular absorption features 
that are readily apparent even at low spectral resolution (R ~ fOOO). For a fixed population 
age, the strengths of the absorption features will depend not only on the overall metallicity, 
Z, but also on the detailed elemental abundance pattern. In general, this greatly complicates 
interpretation of absorption features and spectral indices. However, with the aide of models 
that allow for variation in both metallicity and abundance pattern, one can search for 
combinations of features that are rel atively robust against abundance patter n variations. 
Thomas. Maraston fc Bender 1 2003b ) refined a proposal by Gonzalez (1993J) of an index 
comprised of magnesium and iron lines that is almost entirely insensitive to the level of 
a enhancement, [a/Fe], and is therefore a relatively robust tracer of Z. The index is 
defined as [MgFe]'= A /Mg b (0.72 Fe5270 + 0.28 Fe 5335), where Mgb, Fe5270, and Fe5335 
are Lick indices defined in IWorthcv et al.l fl99J). The power of this index, and others 



like it (see e.g., iBruzual fc Chariot 2003), is that variable abundance pattern models are 
then not required to measure total metallicities from galaxy spectra. However, a significant 
caveat to this approach is the assumption that abundance patterns in galaxies can be 
described by two parameters: [Fe/H], and [a/Fe]. In reality, analyses of high resolution 
stellar spect ra in the Galaxy consistently demonstrate that the a elements do not var y 
in lock-step jEdvardsson et al.ll 19931 . 1 Venn et all2004l , lFul~bright, Mc William fc Richll2007r i. 
This means that the true total metallicity, Z, can only be accurately estimated if the detailed 
a bundance patterns ar e known. 



Gallazzi et alj (|2005l ) applied the idea that [MgFe]' is insensitive to [a/Fe] (assuming that 



the a elements, C and N all vary in lock-step ) to th e analysis of 175,000 galaxies from the 
SDSS with the aide of the Bruzual fc Chariot 120031 ) SPS models. Several age-sensitive and 
Z— sensitive spectral indices were used in their analysis. They derived light-weighted stellar 
metallicities for individual galaxies with statistical errors as low as 0.1 dex for the highest 
S/N spectra. These authors also explored a variety of systematic effects and concluded that 



38 



these could add ~ 0.1 dex systematic uncertainty to the error budget on Z. Probably the 
dominant systematic in the results of Gallazzi et al., and others based on similar techniques, 
is the assumption that the a elements and C and N all vary in lock-step. 

Another technique is to fit the entire optical spectrum with models that allow for variation 
in age and total metallicity, and to simply ignore the possibility that the observed spectrum 
may have non-solar abundance patterns. This appr oach has become incr easingly common as 
automated spectral fitting codes garner wider use. iToieiro et al.l (|201lT l fit stacked spectra 
of luminous red galaxies from SDSS as a function of luminosity and redshift with the VESPA 
code. Remarkably, they found essentially no variation in Z with g alaxy luminosity (though 
the sa mple was restricted t o very luminous galaxies), and for the Conrov. Gunn fe White! 
( 20091 ) and lMarastonl |2005l ) SPS models there was no evolution in Z with redshift over the 
interval 0.1 < z < 0.5. When the lBruzual fc Chariot] (|2003l ) models are employed, VESPA 
favored higher metallicities at higher redshift. Variation in the best-fit metallicity between 
S PS models po i nts to ward systematic uncertainties in the models at the ~ 0.2 dex level. 

Koleva et al. ( 20081 ) investigated the accuracy of full spectrum fitting codes at recovering 



metallicities. They considered both the ability of one SPS model to recover the known 
metallicity of a second model and the ability of an SPS model to recover the known metal- 
licities of Galactic globula r clusters. They concluded that the high-resolution version of the 



Bruzual fc Chariot |2003h model carries large systematic uncertainties due to its use of the 



STEL1B stellar library. STELIB does not sample a wide enough range in metallicity, and so 
metallicities derived from the Bruzual & Chariot models will carry large uncertainties. The 
other models tested, including the Vazdekis/MILES and Pegase-HR models, reproduced 
the known globular cluster metallicities to an accuracy of « 0.14 dex. 

Of course, the most robust method for estimating stellar metallicities is with models that 
self-consistently include abundance pattern variation. Such models have existed for the 
Lick indices since the late 1990s while models for the full optical-NIR spectrum have only 
recent ly become available (jCoelho et al. 2007 . Walcher et al. 20091 . Conrov fc van Dokkum 



2012a). 



Trager et al 



2000bl laD modeled three spectral indices with SPS models that allowed for 
variation in abundance patterns. Trager and colleagues focused on fitting for the age, total 
metallicity, Z, and enhancement ratio, [a/Fe]. In practice, they grouped elements into two 
types, the 'depressed' and the 'enhanced' group, with the latter comprised of a elements 
except Ca, plus C and N, while the f ormer comprised the iron-peak elements and Ca. In 
subsequent work. lThomas et al. modeled the same set of spectral indices as in Trager 

et al., though with their own SPS models and with an expanded dataset. Their conclusions 
largely echoed those of Trager and collaborators, with the novel finding that at fixed velocity 
dispersion galaxies in denser environments appeared to be more metal-rich than galaxies in 
less dense environments. 

As discussed i n Sec tion [ 4] most analyses of Lick ind ices fit single-age models to the data. 



Serra fc Trager (2007) and TTrager fc Somerville ( 20091 ) considered the effect of multiple age 



components in the fitting of the Lick indices of mock galaxies and found that the derived 
metallicities agreed well with the light-weighted mean metallicities when the youngest com- 
ponent was older than ~ 3 Gyr. For young components with ages 1 — 3 Gyr there was 
a slight bias in the sense that the derived metallicities were ~ 0.1 dex higher than the 
true light-weighted mean metallicities. These authors did not consider young components 
younger than 1 Gyr. 

The estimation of total metallicity with models that include abundance pattern variation 
are the most robust but also the most restricted. Models such as those by Thomas et al., 
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Conroy & van Dokkum, Worthey, or Schiavon are limited to ages > 1 Gyr, and therefore 
are routinely applied only to early-type galaxies. A major area for future growth in this 
field is the extension of these models to younger ages in order to model galaxies of all types. 



5.4 Abundance Patterns 

Wallersteinl (| 19621 ) was the first to show that non-solar abundance patterns were common 
amongst metal-poor stars in the Galaxy. It was therefore natural to wonder if other galaxies 
harbored non-solar abundance patterns. Early d ata suggested that magnesium may vary 
more than iron amongst the early- type galaxies I O'Connel] 19761 . Peterson! 19761 . Peletier 
19891 ). b ut it was with the availability of new , relatively high-quality spectral synthesis mod- 
els that IWorthev. Faber fc Gonzalea (|l992T l were able to firmly establish that the [Mg/Fe] 
ratio varied within the early-type galaxy population. The discovery of non-solar abundance 
patterns has been one of the most significant discoveries afforded by population synthesis 
models. The ratio between the abundance of a— elements and iron-peak elements, [a/Fe], 
is considered to be particularly valuable, as it provides insight into the SFH of galaxies 
( Tinslevl[l979l ). following the development of models that allowed for arbitrary variation 
in abundance patterns, the measurement of abundance ratios became routine. Particular 
attention has been paid to [Mg/Fe], in large part because it is the most readily measurable 
abundance ratio in moderate resolution spectra, thanks to the numerous strong Fe I features 
and the Mgb feature a t 5200 A (see Figure [ 3 



As discussed above, iTrager et al 



(2000b a) presented the first systematic investigation 
of abundance patterns in early-type galaxies. They found enhancement ratios, [a/Fe] gen- 
erally in excess of zero, with a strong correlation with velocity dispersion. As Trager and 
collaborators only fit H/3, Mgb, and two iron indices, their [a/Fe] is in reality a statement 
about [Mg/Fe]. Since then, numerous authors, using different SPS models, have all reached 
the s ame conclusion that [Mg/Fe] exceeds zero and increases with galaxy velocity disper- 



sion llThomas et al.|[20oi.lSchiavonl2007l . ISmith et al1l2009l . iGraves. Faber fc Schiavonll2009l . 



Johansson. Thomas fc Maraston 2012h . 

As noted in the previous section, early studies grouped elements into two classes, 
hanced' and 'depressed', and so only the total me t allicity and the enha ncement factor could 



2000al . iThomas et al.ll2005l ). More recent work 
has allowed for the abunda nce of a larger number of elements to vary, including C, N, Ca, 



be measured from the data (jTrager et al 



and Ti (see Schiavon 2010l . for a complementary review to the one presented below). 

The calcium abundance is most commonly measured from the Ca I line at 4227 A, as 
probed by the Ca4227 Lick index (see Figure |3}- Other regions of the spectrum with 
strong sensitivity to calcium abundance include the Ca II H&K lines and the triplet of 
Ca II lines at 8600 A (CaT). The calcium abundance in early-type galaxies has been de- 
scribed as a 'puzzle' because its abundance does not seem to track magnesium, despite 
th e fact that they are b oth a elements. This w as noted qualitatively in the early work 
of IVazdekis et alj ( 1997 ) and Trager et al. I dl998h and quantified with variable abundance 



models by Thomas. Maraston fc Benderl I 2003ah . The CaT feature is weaker than even 
solar-scaled models, for reasons t hat w ill be discussed in Section [7] Schiavon] ( 2007 ) and 
Johansson. Thomas fc Maraston ( 20121 ) fit the Ca4227 index (along with others) of early- 



type galaxies from SDSS and find [Ca/Fe] values close to zero. 

Carbon and nitrogen abundances can be measured from the NH, CH, C2, and CN 



molec ular absorption features in the blue spectral region (see Figure [3]). Kelson et al 



( 20061 ) were the first to probe the carbon and nitrogen content of early-type galaxies 
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based on Lick indices. ISchiavonl (|2007l ) measured [C/Fe] and [N/Fe] from Lick indices as 
a function of early-type galaxy luminosity based on stacked spectra from SDSS. He found 
that the abundance of these elements increased significantly with increasing luminosity. 
Johansson. Thomas fc Maraston reached similar conclusions based on analysis of in- 

dividual early- type galaxies in SDSS. The NH feature at ~ 3360 A is in principle a powerful 
probe of [N/Fe] because it is sensitive only to N (unlike CN). To date this feature has not 
been extensively utilized, both because it is near the atmospheric transmission cutoff and 
thus difficult to observe and because variable abundance models are not well-developed in 
this spectral region. Nonetheless, several authors have investigated empir ical trends of this 
feature with other properties of early- type galaxies, with mixed results (jDavidge fe Clark 



1994. Ponder et al 



properties or early-type f 
ll99Sl . lServen et al.ll201lh 



Moderate resolution spectra of early-type galaxies respond strongly to changes in the 
abundances of Fe, Mg, C, N, Ca, and Ti, and thus these elements are the most easily 
measured. However, the optical and NIR spectrum of an old stellar population contains 
i nform ation on many more elements. In a forward-looking study. IServen. Worthev fe Brilev 
(2005) considered the spectral response due to variation in the abundance of 23 separate 
elements and concluded that for S/N= 100 spectra and a galaxy velocity dispersion of 
o = 200 kms -1 , one could in principle measure the abundance of C, N, O, Na, Mg, Al, Si, 
Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Sr, and Ba. Of course, as the effects are often quite subtle, 
all relevant systemat ics must be well-controlled. As a firs t step toward extracting these 
more subtle features, IConrov. van Dokkum fc Graves! (|2013h reported measurements of the 
abundances of the neutron-capture elements Sr and Ba from analysis of moderate-resolution 
spectra of nearby early-type galaxies. Studies such as these should provide important new 
clues to the formation histories of old stellar systems. 

Attention here has focused exclusively on early-type galaxies with little or no ongoing 
star formation, and for good reason. The hot massive stars that dominate the SEDs of 
actively star-forming galaxies have generally weak metal absorption lines, especially in the 
optical-NIR, and so their presence acts primarily to dilute the strength of the metal lines of 
the older, cooler stars. More fundamentally, models simply do not currently exist to study 
the abundance patterns of actively star-forming galaxies. 



What About Oxygen? Oxygen is the most abundant element in the universe after 
hydrogen and helium. Its abun dance relative to iron, [O/ Fe] , is therefore a critical variable 
in both stellar interior models (jVandenBerg fc Bell 2001) and stellar atmospheres. There 
are no transitions of atomic oxygen visible in moderate resolution spectra, but molecules 
involving oxygen (including TiO, H2O, CO, and OH) do create significant features in the 
optical-NIR spectra of cool stars. TiO is the only molecule with transitions in the op- 
tical (beginning at ~ 4500 A), and unfortunately the work on response functions of the 
Lick indices to abund a nce changes have neglected TiO l ine lists from their computations 
( Tripicco fc Bell! Il995l . iKorn. Maraston fc Thomas! 120051 ). The situation with oxygen is 



quite complex because CO has the highest dissociation potential of any molecule, and thus 
increasing the abundance of oxygen has a cascading effect on other species. For example, in- 
creasing oxygen causes more carbon to be consumed by the formation of CO, which in turn 
causes a decrease in the concentration of CN and C2. The oxygen abundance must therefore 
be known, or assumed, before one can interpret other spectral features. The OH lines in the 
ff-band probably afford the most direct constraint on the oxygen abundance, but these fea- 
tures have not been exploited because of a lack of available models. This limitation has now 
been overcome with the variable abundance ratio models of lConrov fc van Dokkuml l|2012ah 



41 



that extend to the Jt-band. The upshot is that direct constraints on the photospheric abun- 
dance of oxygen in other galaxies currently do not exist, but they are possible with newly 
available models. This presents a significant limitation to current models because of the 
strong effect of oxygen on the effective temperatures, etc., of stars. More troubling is the 
observation that in the Galactic bulge the abundan ces of oxygen and magnesium do not 
track each other (jFulbright. Mc William fc RichlfeoOTT l. In fact, bulge stars with [Fe/H]> 0, 
have [0/Fe]< 0.2, [Si/Fe]< 0.2, [Ca/Fe]< 0.2, and [Ti/Fe]< 0.2, while [Mg/Fe]~ 0.3 for 
the same stars. The lesson to draw from the Galactic bulge is that one should be very 
cautiou s in using [Mg/Fe] as a tracer of [a/Fe]. Because of these uncertainties surrounding 



oxygen, 



Thomas et al.l (|2005f l investigated its effect on their analysis of Lick indices. They 
concluded that if oxygen is in reality decoupled from the other a— elements (e.g., [O/Mg] 
is non-zero) then the ages, [Fe/H], and [Mg/Fe] abundances are largely unchanged with 
respect to their standard model where oxygen tracks the other a— elements. However, these 
authors only considered the effect of oxygen on the stellar spectra: they did not consider the 
non-trivial effect of oxygen on the isochrones (jDotter et alj 120071 ). Further work is clearly 
needed regarding the important role of oxygen in SPS models. 



Z vs. [Fe/H] The iron abundance, [Fe/H], is often regarded as a proxy for the total 
metallicity, Z, but, as discussed above, the correspondence be tween these two quan tities 



breaks down for non-solar abundance patterns. As shown in iTraeer et all (|2000bh . the 
relation between [Z/H] and [Fe/H] can be parameterized by [Z/H]= [Fe/H] [a/Fe]. The 
constant A depends on the detailed abundance pattern; for models in which all a elements 
and C and N are grouped together, A — 0.93, while for models with [0/Fe] = [C/Fe]= 0, 
A = 0.77. It needs to be emphasized here that formulae such as the one presented above 
only apply in the simple case in which all a elements vary in lock-step and C and N either 
are solar-scaled or vary with a. This is probably not a good approximation to reality. From 
the mo deling perspective, some have chosen to enhance or depre ss element groups at fixed 
Z (e.g. jTrager et al.ll2000bl iThomas. Maraston fc B ender 2003b), while others have chose n 



to fix [Fe/H] as other elements vary (e.g., Schiavonf 2007 . Conrov fc van Dokkum 2012a ) 
The choice would be rather arbitrary and inconsequential were it not for the fact that 
most models adopt a fixed set of isochrones as the abundances are varied. The important 
question is then wh ether isochro n es va ry more at fixed Z or fixed [Fe/H]. At least for 
variation in [O/Fe], iDotter et al.l (|2007l ) have demonstrated that the isochrones vary in 



T e g by approximately the same absolute amount (though with a different sign) when the 
variation is computed at fixed Z or fixed [Fe/H]. The disadvantage, at least conceptually, of 
computing models at fixed Z is that a variation in one element forces variation in another 
so that Z is held constant. 



5.5 The Low-Metallicity Tail 

Stars within the Galaxy have a wide range of metallicities, and so one might expect other 
massive galaxies to harbor stars of varying metallicity. This expectation was confirmed 



with resolved color magnitude diagrams of nearby early- type galaxies (e.g.. iGrillmair et al 



19961 . lHarris fc Harrisll20o"ol iMonachesi et al.ll201ll ). The difference between a monometallic 



metal-rich stellar population and one that includes even a small fraction of metal-poor 
stars will be most apparent in the UV. This occurs for two reasons: at low metallicity the 
main sequence turn-off point is considerably hotter than at solar metallicity, and at low 
metallicity stellar populations are observed to contain significant numbers of blue HB stars. 
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The closed box model of galactic chemical evolution predicts that ~ 10% of the stars in 
an evolved galaxy will have [Fe/H]< — 1. Old stellar populations with [Fe/H]< —1 can be 
expected to harbor blue HB stars (by analogy with the Galactic globular clusters), and 
thus there will be a large difference in the UV between a closed box metallicity distribution 
fu nction and one that is monometa llic . 

Bressan. Chiosi fc Fagotto! (1994) constructed spectrophotometric evolutionary models 



for elliptical galaxies and found that closed box models produced an excess of UV flux com- 
pared to observations. They argued that the failure of the model was due to the excess of 
metal-poor stars and concluded that real early-type galaxies contain a more narrowly peaked 
metallicity distribution function. This conclusion mirrors the observation in the solar neigh- 
borhood, where fewer low-metallicity stars are observed than predicted by the closed box 
model (this is known as the 'G-dwarf problem'; Tinslev 1980h . Worthev. Dorman fc Jones] 
(1996) arrived at a similar conclusion from analysis of a UV-optical spectrum of M31 and 
spectral indices for a sample of early-type and SO galaxies. These authors argued for a 
metal-poor component a factor of 2 — 3 lower than predicted by th e closed box m odel. 
From the analysis of Lick indices of a sample of early-type galaxies, iGreggid (1 19971 ) also 
concluded that a closed box model for the metallicity distribution failed to reproduce the 
data. 

While it is reasonably clear that early-type galaxies contain fewer low-metallicity stars 
than predicted by the closed box model, the more general task of measuring the detailed 
metallicity distribution function, or even the fraction of low-metallicity stars from integrated 
light measurements is much more complex. As discussed in Section 14.51 the UV spectral 
region is sensitive to minority populations of hot stars of many types, including metal-poor 
main sequence stars, young stars, and exotic and uncertain evolutionary phases such as blue 
stragglers, extreme HB stars, and post-AGB stars. The latter two are very hot, and so are 
more influential at A < 2000 A, where the UV upturn phenomenon in early- type galaxies 
exists. From co nsideration of the mid-UV (2000 A < A <4000 A) spectra of M31 and other 
nearby galaxies, Worthe y, Dorman fc Jones! ( 19961 ) argued for metal-poor fractions < 5%. 
Maraston fc Thomas! (|2000h mo deled the mid-UV spectra of four nearby galaxies and found 



metal-poor fractions < 6%, and Lotz. Ferguson fc Bohlir] J2OO0I') argued in favor of small, 
but non-zero fractions of metal-poor stars. Maraston fc Thomas (|2000l ) also highlight the 
point, first raised by IGreggid (|l997l ). that known metallicity gradients within galaxies, 
when combined with projection ef fects, will 'contaminate' nuc lear spectra of galaxies with 
metal-poor stars. It is notable that Maraston fc Thomas ( 2000l ) find evidence for a spread in 
metal-poor fractions, ranging from — 6%. This suggests that one should include the metal- 
poor fraction as an additional free parameter when fitting the mid-UV of early-type galaxies. 
Care must be taken when interpreting these fractions because the definition of 'metal-poor' 
varies, although it generally refers to metallicities at least as low as [Fe/H]< —1.5. The 
problem boils down to the fact that the hot star fraction can be measured from UV (and 
perhaps blue spectral) data, but it is very challenging to convert this hot star population 
into a fraction of metal-poor stars owing to uncertainties in advanced phases of stellar 
evolution. 

For actively star-forming galaxies essentially nothing is known regarding their internal 
distribution of metallicities, as the young stars overwhelm the subtle UV features signaling 
the presence or absence of low-metallicity stars. 
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Figure 12: 

The IRX-/3 relation at z = from the Local Volume Legacy Survey (LVL; iDale et aljliooil . left 
panel) and at z ~ 2 from lReddv et al] ll201dl , right panel). The starburst (Meurer) relation is 
shown as a dotted (solid) line in the left (right) panel. In the left panel, symbols indicate 
starburst and star-forming galaxies, along with galaxies split by morphological type. In the right 
panel, symbols indicate galaxies selected in several different ways, either as having a young 
light- weighted age, high bolomctric luminosity, or on the basis of x-ray or 24 /im fluxes. There is 
an overall trend that redder UV slopes correspond to higher IRX values. However, the more 
striking impression is the tremendous scatter in IRX at fixed /3, especially in the region of the 
diagram where most galaxies reside, at —2 < /3 < — 1. 

5.6 Summary 

The technique of measuring metallicities and abundance patterns from old stellar systems 
is well-established. The age-metallicity degeneracy can be robustly broken by considering 
the hydrogen Balmer lines in conjunction with iron features in moderate resolution optical 
spectra. Abundance ratios including [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] are now rou- 
tinely measured from optical spectra. Other elements including O, Na, Si, Sc, Ti, V, Cr, 
Mn, Co, Ni, Sr, and Ba should also be measurable from high S/N optical-NIR spectra of old 
stellar systems, and models are currently being developed to extract this information. It is 
possible to measure metallicities based on photometric data, but the derived values depend 
on the underlying SPS model and on various assumptions like the form of the SFH. Photo- 
metric metallicities are therefore probably accurate only in a relative sense. Spectroscopic 
metallicities seem to be accurate at the ~ 0.2 dex level, with model systematics again being 
the dominant source of error. There is clearly information on the distribution of stellar 
metallicities within galaxies, and previous work suggests that it is possible to measure at 
least the low metallicity tail of the distribution. Further modeling is required to understand 
precisely how much information on the metallicity distribution function is available in the 
SEDs of quiescent galaxies. 



44 



6 DUST 



6.1 Total Dust Attenuation 

There exist four primary techniques for measuring dust attenuation in galaxies. (1) As the 
obscuring effects of dust increase with decreasing wavelength, multi- wavelength analysis 
(including moderate-resolution spectra) of the UV-NIR SED can yield constraints on the 
total dust attenuation. (2) One can exploit the fact that the Balmer line ratio Hq/H/3 
(sometimes referred to as the Balmer decrement) in the absence of dust can b e computed 
from first principles and depends only weakly on the temperature of the gas 1 Qsterbrockl 



198 



^). Comparison of the observed to the intrinsic ratio provides a measurement of the 



attenuation toward H II regions. (3) Energy conservation demands that whatever photons 
are absorbed in the UV/optical must be re-radiated in the IR. The ratio of IR to UV lumi- 
nosity (often referred to as 'IRX') therefore provides a measure of the total dust absorption. 
(4) A sufficiently luminous background source with a known intrinsic spectrum can pro- 
vide a constraint on the attenuation in the foreground galaxy. This last technique will not 
be discussed herein, as it falls outside the scope of what one can learn from the SEDs of 
galaxies. 

Techniques (l)-(3) provide complementar y constraints on the dust within galaxies. Th is 



was most clearly demonstrated in the work of lCalzetti. Kinney &: Storchi-Bergmannl (| 19941 ) 



who analyzed the reddening in a sample of 39 starburst and blue compact galaxies from 
both the Balmer line ratio Hq/H/3 and the slope of the UV continuum. Calzetti and 
collaborators found that the dust optical depth measured from the Balmer lines was nearly 
twice that estimated from the continuum (sometimes referred to as the '1:2' rule). They 
speculated that this could occur if young hot stars are embedded in d usty natal cloud s while 
older stars are subject to attenuation only from the diffuse ISM. Chariot fc Falll ( 200fj| ) 



later provided a detailed model for dust absorption which adopted this physical picture 
wherein young stars are embedded in their natal dusty clouds and, in addition, all stars 
experience dust attenuation due to the diffuse ISM. In this two-component model one 
must specify a transition time, after which the natal cloud is dispersed (typically taken to 
be 10 7 yr), and attenuation curves for both the natal cloud and the diffuse ISM. In the 
standard implementation, both curves are taken to be power-laws with the same exponent 
(jd cx A~ 0,7 ) and the normalization of the natal cloud attenuation curve is taken to be twice 
that of the diffuse ISM. The Chariot & Fall dust model is widely employed when modeling 
dusty SEDs. 

In subsequent work, iMeurer. Heckman fc Calzetti 1 19991 ) compared the UV cont nuum 



slope, ft, to the ratio between FIR and UV flux, log(Lm/ Luv) =IRX, and concluded that 
the starburst galaxies formed a well-defined sequence in the IRX-/3 plane (see Figure [12]). 
With additional assumptions about the underlying stellar populations and dust attenuation 
curve, these authors were able to relate IRX (and hence /3) to the total FUV dust atten- 
uation. The resulting 'Meurer relation' has experienced widespread use for estimating the 
dust attenuation in distant galaxies for which only measurements of f3 are available. How- 
ever, many subse quent studies hav e cast doubt on the practice of using f3 to estimate UV 
dust attenuation. iBell et aL 1 2002h noted that the relationship between /3 and UV attenu- 
ation measured in H n regions in the LMC differed from the Meurer relation. iKong et al 



(2004) demonstrated that, while dust content was the main variable responsible for the 
Meurer relation, other stellar population parameters, particularly the recen t SFH, drive 
galaxies to other locations in the IRX-/3 plane (see also Boauien et al. 20121 ). The under- 



lying extinction curve and star-dust geometry also play an important role in determining 
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where galaxies lie in this plane ( Chariot fe Falll 2000l . Kong et al. 20041 . Panuzzo et al. 2007 , 
Conrov. Schiminovich fc BlantorIl2*01C ) . 



Based on UV-IR photometry of 1000 galaxies. I Johnson et al.l (|2007l ) cast further doubt on 
the utility of the IRX-/3 relation for estimating dust content. Their conclusion rested partly 
on the observation that the relation defined by normal star-forming galaxies is so steep that 
even small errors in /3 will translate into large errors on I RX. Figu re 1121 shows the IRX-/3 
relation for galaxies from the Loca l Volume Leg a cy Su rvey jPale et aliboosh . and the IRX- 
13 relation for z ~ 2 galaxies from iReddv et al. |2oio| y It is clear that, while galaxies of a 
particular type define their own IRX-/? relations, the overall galaxy population presents a 
complicated picture. Fortunately, as multi-wavelength datasets become increasingly com- 
mon, the estimation of dust content from crude measures such as the IRX-/3 relation will 
give way to more robust techniques. 

SED fitting of UV-NIR data has become remarkably common in the past decade. In 
nearly all cases the dust content is a parameter included in the fit, although the derived 
dust opacities are rarely discussed. The reason for this seems to be the fact that the dust 
opacity is a relatively poorly constrained quantity when data only reach the NIR. As an 



examp le of the difficulty in constraining the dust opacity from broadband data, ISalim et al 



(2007J) fit UV-NIR photometry to models and were able to recover V— band dust opacities 
with an uncertainty of ~ 0.4 mag. Furthermore, there was only a weak correlation between 
the dust opacity deriv e d from broadband data and the opacity derived from emission lines 



Brinchmann et al 



1 20041 ) . Some of the difference can be attributed to aperture effects 
and real differences in the opacity measured from emission lines and stellar continua, but 
overall the weak correlation confirms that dust content measured from broadband data 
is quite uncertain. Poorly constrained dust opac ities appears to be a generic result of 
the modeling of UV-NIR broadband SEDs (e.g. 



Shaplev et alll2006l ; iKriek et al]|2008l ; iTavlor et al 



Papovich. Dickinson fc Ferguson! 2001 



20111 ). 



Broadband UV-NIR data do not provide strong constr aints on dust attenuation because 
of th e well-known degeneracy between age and dust (e.g.. iPapovich. Dickinson fc Ferguson 
2001). This degeneracy can be broken, and stronger constraints on the dust opacity can be 



obtained, when either spectroscopic or FIR data are available. Spectroscopic data offer an 
estimate of the mean stellar age and metallicity independent of dust content because narrow 
spectroscopic features a re relatively immune to overall changes in the continuum shape. 



Kauffmann et alj |2003l ) exploited this technique to measure mean ages of ~ 10 5 SDSS 



galaxies from H5 and the 4000 A break. With the ages thus constrained, an estimate of the 
dust attenuation could be obtained by comparing the model colors to the observed ones. 
This study was limited primarily by the fact that the SDSS spectra sample a modest fraction 
of the total galaxy, and so a correction must be made for this aperture effect (ongoing and 
future planned IFU surveys such as CALIFA and MaNGA will remedy this limitation). 
Narrowband photometric surveys also offer the promise of separately constraining age and 
dust attenuation, primarily because of the strongly age-sensitive _D n 4000 spectral feature , 
which is so strong that it ca n be probed with low resolution data (e.g., Whitaker et aLlfeoid . 



Perez-Gonzalez et al. 20121 ). FIR data can also provide a robust measurement of the total 
dust attenuation because they significantly reduce the degenera cy between age and dust 



(e .g.. iBurgarel 



Wuvts et al 



.. Buat fc I elcsias-P aramdl200 
(|2007l ) and'" 



y rec 
20051: 



Noll et al 



2009a). 



Williams et all (200 



NIR color-color diagrams (specifically U — V vs 

diagram) to efficiently separate quiescent and star-forming galaxies 
has since become a popular tool for separating these two galaxy types, especially at high 



ij) pioneered the use of restframe optical- 
V — J; commonly referred to as the UVJ 
The UVJ diagram 
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redshift where data are limited. In a broad sense, the UVJ diagram allows one to break the 
age-dust degeneracy, but only in the sense of separating quiescent and dusty star-forming 
galaxies. Within the star-forming sequence, the age-dust degeneracy is still manifest when 
only broadband information is available. 

Extre mely blue UV c o ntinu a can place strong cons t raints on the amount of dust in the 
system. Bouwens et al. ( 20121 ) and Finkelstein et al. (2012) measured the UV continuum 
slopes of high-redshift galaxies and found that ft systematically decreases toward higher 
redshift and fainter galaxy luminosities. At the highest redshifts and faintest luminosities 
they find —2.0 < ft < —2.5, which not only suggests that these galaxies are dust-free but 
also begins to place interesting constraints on their stellar populations. The models of 



Schaerer! (|2003h reach such blue slopes only at very young ages (~ 10 6 yr) and/or very 



low metallicities (< 10~ Zq), and the inclusion of nebular continuum emission, which 
presumably must be present whenever stars of these young ages are still alive, sets a lower 
limit to the UV continuum slope of —2.5 < ft < —3.0, depending on the details of the SFH. 
Clearly more data probing the restframe UV at very high redshifts would be valuable. 



6.2 Constraints on the Attenuation Curve 



Calzetti. Kinney fc Storchi-Bergmannl (|l994T ) estimated a mean attenuation curve from 39 



starburst and blue compact galaxies based on UV-optical spectra. By assuming that their 
sample was comprised of galaxies with the same underlying stellar populations, they were 
able to construct average attenuation curves by comparing galaxies with high and low 
Balmer line ratios, i.e., by comparing dusty to dust-free galaxies. The resulting average 
attenuation curve, measured over « 1200 A-8000A, was greyer (shallower) than either the 
MW or LMC extinction curve and lacked the broad 2175 A dust feature visible in the MW 
and LMC. The derived curve is now referred to as the "Calzetti attenuation law" and is 
widely used when modeling the broadband SEDs of galaxies. 

The Calzetti law was estimated from a relatively small sample of unusual (i.e., starburst 
and blue compact) galaxies at 2 = and so one may wonder to what extent it is applicable 
to other systems. As noted in the previous section, the IRX— ft relation is sensitive to the 
attenuation curve. The large range in ft at a fixed IRX, evident in Figure 1121 suggests 
considerable variation in the attenuation curve, but interpretation of the IRX— ft plane is 
complicated by the additional variables affecting IRX and ft. 

Progress can be made by estimating stellar population parameters such as the SFH and 
metallicity independently of the dust content (e.g., via narrow spectral features). With 
such quantities in hand, one can exploit the technique pioneered by Calzetti and collabora- 
tors by computing r atios of SEDs of more to less dusty galaxies with the same underlying 



stellar populations. iJohnson et alJ (|2007l ) estimated mean attenuation curves in this way 
for galaxies as a function of the 4000 A break strength and stellar mass. The curves were 
constrained from broadband FUV-NIR photometry and were consisten t both with a simple 



power-law, tx A , and with the Calzetti law. IWild et al.l (|201ll ) exploited a similar 
technique using pairs of galaxies with similar gas-phase metallicities, specific star formation 
rates, and inclinations, and significantly different Balmer decrements to construct attenua- 
tion curves from broadband FUV-NIR data. These authors found significant variation in the 
mean attenuation curve as a function of specific SFR, inclination, and stellar mass surface 
density, based on a sample of 20,000 2 ~ star-forming galaxies. They also found variation 
in the ratio between dust opac ity in the stellar continuum and in the nebula r line emission, 



implying that the '1/2' rule of Calzetti. Kinnev fc Storchi-Bergmannl 1 19941 ) — one unit of 
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continuum opacity for ever two units of line opacity — is not universal. Buat et al, ( 20121 ) 
modeled the UV-FIR SEDs of galaxies at 1 < z < 2, allowing for variation in the dust 
attenuation curve. They found evidence for a steeper attenuation curve (i.e., faster rise in 
the UV) than the Calzetti law in 20 — 40% of their sample. 

6.2.1 The 2175 A dust feature The absence of the 2175 A dust feature from the 



Calzet ti attenuation law is striking given its ubiquity in the MW and LMC. IWitt fc Gordon 
( 20001 ) were able to reproduce the Calzetti law in their radiative transfer models only by 
assuming SMC-type dust, i.e., by ado pting an unde r lying extinction curve that lacked 
the 2175 A dust feature. In contrast, 



Granato et al 



(|2000l ) were able to reproduce the 
Calzetti law in their radiative transfer models with MW-type dust. In their model the 
star-dust geometry plays a central role in shaping the starburst attenuation curve. Their 
molecular clouds are optically thick and so the attenuation curve is deter mined primarily 



by th e fraction of stars ins i de mo lecular clouds as a function of age (see also lPanuzzo et al 



20071 ). iFischera fc Dopital |201ll ) were able to reproduce a Calzetti-like attenuation curve 



if the carrier responsible for the 2175 A feature (e.g., PAHs) was destroyed at high column 
density. 

In all radiative transfer dust models the expectation is that normal star-forming galax- 
ies should show evidence for the 2175 A dust feature provided that the underlying grain 
population is similar to the MW or LMC. Owing to a paucity of re stframe UV spectra of star- 
formin g galaxies, it has proven difficult to verify this expectation. Burgarella. Buat fc Iglesias-Paramol 
(120051 ) analyzed the UV-FIR SEDs of 180 star-forming galaxies and concluded that the 
data required an attenuation curve with a 2175 A dust featur e that was on average weaker 
than observed in the MW extinction curve. I Noll et al.l (|2009bl ) presented stacked restframe 
UV spectra of z ~ 2 star-forming galaxies and found st rong evidence for the presence of 
the 21 75 A dust feature in a subsample of their objects. IConrov. Schiminovich fc Blanton 
( 20101 ) analyzed UV-NIR photometry of z ~ star-forming galaxies as a function of in- 
clination and found evidence for a 2175 A dust feature w i th a strength slightly weaker 
than observed in the MW extinction curve. iHoversten et al.l (|201ll ) analyzed medium band 
UV photometry from the Swift UV/Optical Telescope, and found evidence for a promi- 
nent 2175 A bump in M81. IWild et al.l ( 201 In found evidence in their derived attenuation 
curves for the presence of the 2175 A dust feature, wit h a slight tendency for ga laxies with 



higher specific SFRs to have lower bump strengths. IWiiesinghe et al.l l|201ll ) compared 



dust-corrected SFRs estimated from Ha and the UV and found agreement only when the 
2175 A feature was removed from the attenuation curve, but these authors did not con- 



sider variation in the shape of the attenuation curve. iBuat et al.l (|2012l ) found evidence 
that the 2175 A dust feature is present in 20% of their sample, with an amplitude that is 
half the strength found in the MW extinction curve. Similar to Wild et al., they found an 
anti-correlation between bump strength and specific SFR. The emerging consensus appears 
to be that the 2175 A dust feature is present in the typical star-forming galaxy over the 
redshift range < z < 2, with a strength that is generally lower than in the MW extinc- 
tion curve, and that increases with increasing inclination and decreasing specific SFR. The 
absence of the dust feature in the Calzetti law may then simply be a reflection of the fact 
that the galaxies in the Calzetti sample have unusually high specific SFRs. An outstanding 
question is whether these trends are due to varyi ng star-dust geometries or changes in the 
underlying grain population. As emphasized in Conrov. Schiminovich fc Blanton ( 2010l ) 
and elsewhere, these results have consequences for interpreting the IRX— j3 relation, since 
/3 is frequently measured over a wavelength range that includes the 2175 A feature. 
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6.3 Physical Dust Properties 



The IR emission by dust grains provides strong constraints on a variety of physical properties 
of a galaxy including the bolometric luminosity, dust temperature(s), dust mass, and relative 
abundance of PAH molecules (or, more generally, the grain size distribution). A constraint 
on the bolometric luminosity helps to break the age-dust degeneracy in the UV-NIR, as 
discussed in the previous section. The other parameters provide clues to the nature of dust 
in galaxies and its spatial relation to the stars and gas. 

It has long been recognized that constraints on the dust mass and dust temperatures 
required data beyond the peak in the thermal dust emission spectrum at ~ 100 /im. The 
SCUBA bolometer array at the James Clerk Maxwell Telescope, operating at 450 iim and 
850 fim, has been extremely influential in this field, and now the SPIRE instrument on- 
board the Herschel Space Observatory, operating at 250 — 500 /J,m, promises to revolution- 
ize the field by dramatically increasing sample sizes over a wide range in redshifts. With 
data obtained f or the SCUBA Local Universe Galaxy Survey, Dunne et al. ( 2000h and 
Dunne fc Ealesl ( 2001 ') constrained dust masses and temperatures for 104 nearby galax- 



ies. The data favored models that contained both cold and warm dust components (at 
~ 20K and ~ 40K, respectively). They also found that dust masses estimated with single- 
temperature models were on average lower by a factor of ~ 2 compared to two-temperature 
models, and the two-temperature dust masses yielded dust-to-gas ratios in agreement with 
m easurements in the Galaxy. 

Draine et al] (|2007l ) presented a thorough analysis of the IR SEDs of galaxies from the 



SINGS survey of 65 nearby galaxies. Thes e authors ana l yzed data from Spitzer, IRAS, and 
SCUBA with the physical dust models of iDraine fc Lil (|2007l ). The model considers grain 
populations exposed to a variety of starlight intensities and includes both thermal emission 
and single photon heating of dust particles. In addition, the fraction of dust mass contained 
in PAH molecules is a free variable. Draine et al. confirmed the results of previous work 
that the sub-mm data, provided by SCUBA, were essential for deriving dust masses with 
an uncertainty of < 50%. They also derived strong constraints on the fraction of dust mass 
in PAHs, the interstellar radiation field strength in the diffuse ISM, and the flux-weighted 
mean radiation field strength averaged over the galaxy. The well-sampled spectral coverage 
from the mid-IR through the sub-mm was essential in order to separately constrain each of 
t hese components. 

Dale et al. (2012) analyzed FIR and sub-mm photometry for SINGS galaxies from Her- 



schel in conjunction with previously available 2MASS, Spitzer, IR AS, ISO, and SCUBA 
data. These authors modeled the SEDs with the lDraine fc Lil (|2007l ) dust model, both with 
and without the Herschel data. They found modest differences (factors of « 1.6 on average) 
in the best-fit dust masses, PAH mass fractions, and fraction of dust emission arising from 
the diffuse ISM when the FIR and sub-mm data were either included or excluded. They 
also demonstrated that dust masses estimated via simple modified blackbody fits signif- 
icantly underestimated the dust masses (by factors of ~ 3 in the worst cases), with the 
underestimation becoming more severe as bluer wavelengths were included in the fits. The 
Herschel data are of such high quality that simplistic, modified blackbody dust models are 
n o longer capable of providing adequ ate fits to the data. 



da Cunha. Chariot fc Elbaa (|200S ) developed a phenomeno logical model for the UV-FIR 
SE Ds of galaxies that is b ased on the lBruzual fc Chariot] (|2003l ) SPS models for the starlight, 
the I Chariot fc Falll l |200Ch dust attenuation model, and dust emission comprised of multiple 
components including an empirical spectrum for the PAH emission, warm and cold thermal 
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Figure 13: 



Model fit to the FUV through FIR SED of NGC 337, from Ida Cunha. Chariot fc Elbad fo00$ ). 
The black line is the best-fit model. Also shown is the unattcnuatcd model starlight (blue line) 
and the dust emission separated into dust in the diffuse ISM (solid green line) and dust in birth 
clouds (dotted green line). Figure courtesy of E. da Cunha. 



dust emission, and emission from stochastically heated grains (see Figure [To')) . Based on 
tests with mock galaxies, they concluded that a variety of parameters could be constrained 
within the context of their model, including the temperature of the cold dust component, 
the fractional contribution of the cold dust to the total FIR emission, and the total dust 
mass. They analyzed the same SINGS galaxies as in Draine et al. and derived dust masses 
that agreed to within 50%, demonstrating that dust masses can be reliably measured with 
re asonably well-sampled mid-IR, FI R, and sub-mm data 



Galliano. Dwek fc Chanial ( 20081 ) employed the dust model of Zubko. Dwek fc Arendt 



(2004|) to constrain a variety of parameters including the total dust mass, fraction of 
dust mass in PAHs, fraction of dust luminosity originating from photodissociation regions 
(PDRs), and the gas-phase metallicity. These authors measured these and other param- 
eters for 35 galaxies with optical, IR, and sub-mm photometry and mid-IR spectr oscopy. 
Although no direct comparison was made with the results from lDraine et alj (|2007l ). there 
appears to be a quantita t ive di fference in the derived PAH mass fractions. As discussed ex- 
tensively in lDraine fc Lil (|2007l ). the optical properties of the PAH molecules are uncertain 
and have in many cases been tuned to fit extragalactic data. It would therefore be valuable 
to compare the results obtained from dust models that assume different treatments of the 
uncertain PAH properties. 



50 



6.4 Cosmic Evolution of IR SEDs 



Measuring the cosmic evolution of IR SEDs is important not only for constraining models 
of galaxy evolution but also because local IR SED templates are frequently employed to 
interpret high-redshift data. For example, the local templates are often used to convert 
24 fim observations of high-redshift galaxies into SFRs. If IR SEDs at fixed Lboi evolve 
w ith time the n it wo uld considerably complicate the interpretation of higher redshift data. 
Pope et al. (2006) presented evidence that the SEDs of ULIRGs at z ~ 2 peak at 



local ULIRGs, in agreement with lPope et al 



longer wavelengths than local ULIRGs, implying that t he dust tempe r ature in high-redshift 
ULIRGs is on average ~ 5K cooler than local ULIRGs. Muzzin et al. I |2O10l ) analyzed high- 
quality data for two 2 ~ 2 ULIRG s and found t hat th ese galaxies contained colder dust than 

l|2006l). Muzzin et al showe d that the SED 
Elbazl (|200ll ) templates that 
were an order of magnitude less luminous than their 2 ~ 2 galaxies. In other words, the 
2 ~ 2 ULIRGs had SED shapes th at were similar to local LIRGs. Exploiting the power of 
Herschel data, Hwang et al. ( 2O10l ) estimated dust temperatures by fitting modified black- 
bodies to the data, and found evidence for slightl y cooler dust tempe ratures at z ~ 1, by 



shapes of their ULIRGs were well-fit by the 2 = 1 Chary 



2 — 5 K on average, for galaxies with Lm > W 11 Lq. iDaddi et al.l ( 20071 ) compared restframe 
8 /im luminosities of z ~ 2 galaxies to local templates and found evidence for excess mid-IR 
emission for galaxies with Z/8 Mm > lO 11 !/^. Th i s mid -IR excess problem was confirmed by 



Papovich etafl |2007t ) and later bv lElbaz et"ZI <|2010h with Herschel data. It is noteworthy 
that the problem only appears at 2 > 1.5. Bolometrically luminous high-redshift galaxies 
thus have colder dust and more flux in the mid-IR compared to galaxies at the same Lboi 
at 2 = 0. The mid-IR excess problem has important implications for the common practice 
of estimating SFRs for galaxies at 2 > 2 with Spitzer 24 \xm observations (i.e., restframe 
~ 8 /mi). The sign of the excess is such that SFRs estimated by extrapolating observed 
24 fim flux with local templates will tend to overestima te the true SFRs, by factors of several 
in the most extreme cases (e.g.. iPapovich et al.ll2007h . 



Elbaz et al.l (|2011m investigated the origin of the differences between high and low redshift 
IR SEDs with very deep Herschel observations. They argued that it is the distribution of 
SED types that is changing at high luminosity as a function of redshift. At low redshift, 
high luminosity systems (i.e., ULIRGs) are predominantly star-bursting and compact, while 
at higher redshift ULIRGs are mostly the high-luminosity extension of normal star-forming 
galaxies. Local ULIRGs therefore have hotter dust temperatures and a depressed emission 
component from PAHs in the mid-IR compared to distant ULIRGs because the former 
are compact starbursts whi le the latter are normal, extended star-forming galaxies (see 
also iRuiopakarn et al.ll2012l . who reached similar conclusions). The physical origin of this 
result is not yet fully understood, but a plausible explanation is that compact starbursts 
have harder and more intense radiation fields. The radiation field obviously has a direct 
influence on the dust temperature, and it may also modulate the dust mass fraction in 
PAH s and therefore the mid-IR luminosity ( Voit 1992 . Madden et al. 20061 . Draine et al 
2007m . 



6.5 Summary 

The Calzetti attenuation curve works remarkably well at describing the mean attenuation 
properties of star-forming galaxies over most of cosmic history. However, closer scrutiny 
of the data reveals variation in the attenuation curve of galaxies as a function of galaxy 
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properties, as expected on theoretical grounds. In contrast to the Calzetti attenuation curve 
for starburst galaxies, the 2175 A dust feature is apparently common in normal star-forming 
galaxies, with a strength that is weaker than observed in the MW extinction curve. The 
IRX— /3 relation should be used with caution, as there is evidently no universal relation for all 
galaxies. In addition, the relation is so steep over the range in j3 where most galaxies reside 
that estimating dust opacity in this way is highly unreliable. Dust masses can be measured 
to an accuracy of ~ 50% when FIR and sub-mm data are available. Physical dust models, 
such as the model of lDraine fe Lil 1)20071 ). are capable of extracting a variety of parameters 
from the global SEDs of galaxies, including the fraction of dust mass in PAHs, the fraction 
of dust emission due to PDRs vs. the diffuse ISM, and the typical interstellar radiation field 
strength heating the dust grains (or, alternatively, the typical dust temperature). Further 
work appears to be needed to sort out which derived parameters are robust, and which are 
dependent on uncertain components of the model. Below Lbol ~ 10 n L©, the IR SEDs of 
galaxies evolve little over the interval < z < 2. At higher luminosities there is a marked 
shift such that high redshift ULIRGs have more mid-IR flux and a colder dust component 
that peaks at longer wavelengths. A plausible explanation for this trend is that ULIRGs 
at high-redshift are at the massive end of the normal star-forming galaxy sequence, while 
ULIRGs in the local universe are unusual, compact, star-bursting systems. 



7 THE INITIAL MASS FUNCTION 



The stellar IMF is of fundamental importance for m any areas of astrophysics. The form 
of the IMF is well-constrained in the Galactic disk ( Salpeter 19551 . Miller fc Scald 19791 . 



KroupallioOll . Ichabrierll2003h . However, it is not clear whether the IMF has had the same 
form over all of cosmic time and in all environments. The low mass end is of particular 
importance, as ~ 60 — 80 % of the stellar mass density in the Universe is in the form of 
stars with masses < 0.5 Mq. Low mass stars are very faint and so contribute only a few 
percent to the integrated light of an old population. For this reason, when fitting the SEDs 
of galaxies, an IMF is typically assumed, and so the IMF is one of the largest sources of 
systematic uncertainty in studies of extragalactic stellar populations. 

It has been known since at least the 1960s that dwarf-sensitive and, to a lesser degree, 



the integrated light of old stellar populations (e.g., SDinradl 1962. SDinrad & Tavlorl 


1971. 


Whitfordll977i, Frogel et al.ll978. Frosjel, Persson & Cohenll980lJVazdekis et all 19961. 


Schiavon, Barbuv & Sin 


19971. ISchiavon et al. 19971. ISchiavon, Barbuv & Bruzual A.II2OO0I). This idea is limited to 



quiescent systems because the massive stars associated with ongoing star formation would 
dramatically outshine the faint low mass stars. The fundamental difficulty then, as now, 
is the separation of abundance effects from giant-to-dwarf ratio effects. Both effects can 
c hange the strength of the gravity-sensitive lines. 



Spinrad fc Tavlorl l| 19711 ) concluded that the nuclei of the nearby galaxies M31, M32, and 
M81 are dominated by dwarfs, with an implied M/L — 44. These c onclusions wer e based on 
photoelectric scanner observations reaching to 1 \im. Subsequently, Cohen ( 19781 ) measured 
the IMF-sensitive Na I, CaT, FeH and TiO features in the nucleus of M31 and M32 and 
concluded that an IMF comparable t o that of the solar nei ghborhood was favored, along 
with metallicity enhancement in M31. iFaber fc French] (|l980l l also measured the strength of 



the Na I doublet in M31 and M32 and instead argued for a dwarf- rich IMF, along the lines of 
Spinrad fc Tavlorl (|l97lh . with M/Lb = 28. iFrogel et all l|l978l ) analyzed optical-NIR col- 
ors and the CO bandhead at 2.3 fim for 51 early-type galaxies and concluded that dwarf-rich 
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IMFs were inconsistent with the data and favored m odels with M/Ly < 10. With updated 
detector technology and more sophisticated models, Carter. Visvanathan fc Pickles ( 19861 ) 
argued against dwarf-rich IMFs on the basis of the Na I, CaT, TiO and FeH features in 14 
elliptical and lenticular galaxies. Their results imp lied mass-to-light ratios of M/L ~ 6. In a 
series of papers, Hardy, Coutu re and collaborators ( Hardy fc Couture! 19881 . Delisle fc Hardy 
19921 . ICouture fc Hardvlll993h analyzed the IMF features Na I and FeH for of order a dozen 
galaxies and concluded that dwarf-rich IMFs were strongly disfavored. Measurements of 
the FeH band at 1 pan (also known as the Wing-Ford band) were particula rly important, a s 
this feature becomes very strong for IMF exponents steeper than Salpeter (|Whitfordlll977l ). 

Early work aimed at measuring the IMF from integrated light spectra suffered from se- 
rious limitations, including (1) the lack of accurate stellar evolution calculations across 
the main sequence and through advanced evolutionary phases; (2) the use of empirical 
stellar spectra collected from the solar neighborhood, which implied that the effect of non- 
solar abundance ratios could not be adequately assessed; (3) poor NIR detector technology, 
which made it very difficult to measure red spectra at the sub-percent level precision nec- 
essary to measure the low-mass IMF in integrated light. These limitations contributed 
to the wide range of reported IMF exponents and mass-to-light ratios noted above. In 
the past two decades, each of these limitations has been significantly reduced, though 
not eliminated. Stellar evolution calculations have improved dramatically since the 1970s, 
the resp onse of stellar spectra to elemental abundance variations has been studied in de- 



spe 

tail (e.g..lTripicco fc Bel]l995 jKorn. Maraston fc Thomasll2005l , IServen. Worthev fc Brilev 



2005 , Conrov fc van DokkumI 2012a ) , and NIR dete ctor technology ha s stea dily improved 



The debate over IMF variation was renewed by ISaglia et al.l (|2002f ) and ICenarro et al 
(120031 ) . who independently reported an anti-correlation between the strength of the IMF- 



sensitive CaT index and velocity dispersion for elliptical galaxies. Comparison to SPS 
models with solar-scaled abundance ratios required dwarf-rich IMFs, with an exponent of 



a 4 at the highest dispersions (compared to the Salpeter value of 2.35). ISaglia et al 



(|2002h argued against this interpretation, as such a steep IMF would imply M/L » 40 



well i n excess of newly available dynamical constraints. I Worthev. Ingermann fc Serven 
(|201lr ) investigated several calcium-sensitive spectral features, and concluded that a mod- 
est decrease in [Ca/Fe] with velocity dispersion could also explain the observed trends. 



Worthev. Ingermann fc Served (|201lf ) also noted that such dwarf-rich IMFs would yield 
optical-NIR colors (e.g., V — K) much redder than even the reddest observed early-type 
galaxies. 

An IMF differ ent from the Galactic di s k (oft en referred to as a non-universal IMF) was 
also reported by Ivan Dokkum fc Conrovl (|2010h based on the analysis of high S/N spectra 
extending to 1 \im for eight massive galaxies in the Coma and Virgo galaxies. These data 
were of extremely high quality thanks to new detector technology installed on the Low 
Resolution Imaging Spectrometer (LRIS) instrument at the Keck Telescope. The data were 
interpreted with new models limited to solar-scaled abundance patterns. Dwarf-rich IMFs 
were favored, with an IMF exponent in the range 2.3 < x < 3. It is worth pausing here to 
note that the meaning of 'dwarf-rich' has evolved over the years. In the 1970s dwarf-rich 
IMFs indicated mass-to-light ratios exceeding 20, ofte n reaching as high as 40 . Wit h the 
IMF in the solar neighborhood now well-established |Kroupall200l] . Ichabrierll2003h . the 
reference point has shifted. Now, 'dwarf-rich' IMFs are those that contain relatively more 
low-mass stars than found in the solar neighborhood and would indicate, amongst other 
implications, a non-universal IMF. In this sense, th e Salpeter IMF is 'dwarf-rich', as it 
contains more low-mass stars than the lKroupal (|200 lM IMF for the solar neighborhood. 
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Figure 14: 



Comparison between the best-fit model and observed spectrum for NGC 4621, a massive 
early-type galaxy. Both data and model have been continuum normalized. The S/N is indicated 
as a grey shaded band; the data at < 4000 A were not included in the fit but arc included here for 
completeness. The model has 21 free parameters including the abundance of 11 elements and the 
shape of the low-mass IMF. The rms difference between the model and data is 0.8%. IMF-sensitive 
features include the Na I doublet at 0.82 (im, the Ca II triplet at ss 0.86 fj,m, and the FeH bandhead 
(i.e., the Wing-Ford band) at 0.99 fim. The best-fit M/L is ~ 2 x as large as would be expected 
for a MW IMF, suggesting that the IMF is not universal. From lConrov fc van Dokkuml j2012bl V 



Conrov fc van Dokkuml (|2012bh were the first to analyze high-quality absorption line 



spectra of 34 early-type galaxies and the nuclear bulge of M31 with a model that al- 
lowed for both IMF variation and abundance pattern variation ( in addition to variation 



i n age and several nuisance parameters). The model, presented in lConrov fc van Dokkum 



(|2012af ). allows for variation in the abundance pattern of 11 eleme nts, and is based on the 
latest isochrone tables and the MILES and IRTF stellar libraries. Conrov fc van Dokkum 



( 2012bl ) simultaneously fit the full blue and red galaxy spectra to their models, including 
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the classic IMF-sensitive features Na I, CaT, and FeH, as well as IMF-sensitive features in 
the blue spectral region. An example fit to a massive early-type galaxy is shown in Figure 
1141 They argued that the models favor IMFs that become progressively more dwarf-rich 
(i.e., more bottom heavy) for more massive and more a— enhanced galaxies. Over the full 
sample, the mass-to-light ratios predicted by their best-fit IMFs vary by only a factor of 
3 at fixed age and metallicity. The IMF variation advocated by IConrov fe van Dokkum 



( 2012b! ) is thus much subtler than previous claims for dwarf-rich IMFs. The best-fit mass- 
to-light ratios do not v i olate dynamical constraints. Similar conclusions were reached by 



Smith. Lucev fc Carter! (12012 ) based on t he an alysis of 1 /im spectroscopy of 92 galaxies in 



the Coma cluster, and bv lSpiniello et al. I (|2012l ). who analyzed indices sensitive to Na I and 
TiO for a sample of early-type galaxies. 

Meanwhile, there have been several less-direct arguments favoring top-heavy (or bottom- 
light) IMFs in some types of galaxies and in some phases of a galaxy's evolution. iTinslev 
(1980) suggested that a constraint on the IMF near the turn-off mass could be obtained 
by comparing the evolution of the M/L ratio and color of passively evolving galaxies. 
van Dokkum ( 20081 ) applied this idea to the evolution of galaxies on the fundamental plane, 
and found evidence that the IMF is flatter than the MW IMF at ~ 1 Mq in early-type 
galaxies. However, Ivan Dokkum fc Conrov (|2012l ) demonstrated that when galaxies are 



compared at fixed velocity dispersion, as opposed to fixed stellar mass, the evolution in 
the fundamental p lane is consistent with the standard Salpeter slope at ~ 1 M© (see also 



Holden et allboiol ). Numerous authors have demonstrated that the evolution of the cos- 



mic star formation rate and cosmic stellar mass density ar e mutually inconsistent but can 



be brought into agreement by adopt i ng a top-heavy IMF (jHopkins fc BeacomlboodlDavel 



200. 



i|, Wilkins. Trentham fc Hopkins! 2008). However, recent improvement in the model- 



ing of the data, including luminosity-dependent dust corrections, rising SFHs, and revised 
values of the faint end slope of the mass function have reconciled the cosmic star forma- 
tion and mass densities derived with standard IMFs ( Reddv fc St cidcl 200_^, Papovich et al 



201ll . iBehroozi. Wechsler fc Conrovll2012f ). A third argument in favor of top-heavy IMFs 



is the apparent inability of galaxy formation models to reproduce the observed number 
of luminous submillimeter galaxies (|Baugh et al.ll2005l ). Here too recent observations and 
more sophisticated m odels have significantly reduced the tension between models and data 
( Havward et al.ll2012l ). It thus appears that much of the evidence originally favoring a top- 
heavy IMF at certain epochs and under certain conditions can now be explained by more 
mundane effects. 



7.1 Summary 



Extremely dwarf-rich IMFs are now routinely ruled out both by sophisticated SPS models 
and dynamical constraints, and earlier evidence for top-heavy IMFs has not held up to fur- 
ther scrutiny. However, more modest IMF variations appear to be supported by the data, at 
the level of a factor of 2 — 3 in M/L. It is intriguing that constraints on the IMF based solely 
on dynamical modeling also favor a mild steepening ( i.e., relatively mor e low- mass stars 



with increasing velocity d i spersion ( Tmi_et_al 



2010. 



Spiniello et"al] I2OI2I. 



Thomas et al 



2012I . button. Mendel fc SlmardlboiH 



2011 



2010. 



Auger et al 



Dutton et al. 



2012a 



201 



Sonncnfcl d et al.l 120121 . ICappcllar i et al 



Graves fc Faber 



b|). As both dynamical and stel- 



lar population-based techniques suffer from non- negligible but largely orthogonal system- 
atics, the most promising direction for future constraints will come from joint analyses of 
the same systems with SPS and dynamical techniques. 
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8 CONCLUDING REMARKS 



The goal of this review has been to summarize what we have learned about galaxies from 
their panchromatic SEDs using the tools of stellar population synthesis. By way of con- 
cluding, I would like to turn toward the future and highlight areas where additional work is 
needed to make SED modeling both a more firmly quantitative science and also an effective 
engine for new discoveries. In short, and not surprisingly, the future requires better data, 
better models, better comparison of models to data, and a better understanding of what is 
in principle knowable from the analysis of galaxy SEDs. 

A clear theme of this review has been the power of combining broadband data with mod- 
erate resolution spectra. The SDSS has been a revolutionizing force in this regard, as it 
has provided high-quality photometry and optical spectra for over one million objects. The 
great drawback of the SDSS is that it is fiber-based and thus the spectra only sample the 
central 3" (in diameter) of galaxies. This drawback will be alleviated with IFU spectroscopic 



surve ys of nearby galaxies, including t he recently completed SAURON survey | Bacon et al 



2001 '). the ongoing ATLAS3D sur vey |Cap pellari et al. 201lh of 260 early- type galaxies, the 
CALIFA survey of ~ 600 galaxies ( Sanchezet"al*ff2012h . and the proposed MaNGA survey of 
a mass-limited sample of ~ 10, 000 galaxies. Such surveys will provide high-quality spectra 
that are well-matched to the broadband data. Restframe NIR spectra will be a new frontier 
for SPS studies as next generation NIR facilities become operational, including ground-based 
spectrographs (FIRE, FMOS, KMOS, MOSFIRE, FLAMINGOS-2, etc.), and the eventual 
launch of the James Webb Space Telescope. High-quality models are only now being devel- 
oped to interpret such data. In addition to spectra, narrow-band photometry will also be 



a val uable addition to the landscape, as pioneered by t he COMBO-17 su rvey (|Wolf et al 



20041 ) and now utilized by surveys such as ALH AMBRA J Moles et al.ll2008h, the NEWFIRM 



Medium-Band Survey (|Whitaker et ailboill ). and SHARDS (|Perez-Gonzalez et all 120121 ). 
Grism data will also help bridge the gap between broadba nd photometry and m oderate 
resolution spectra, as demonstrated by the 3D-HST survey ijBrammer et al.ll2012h . In ad- 
dition to these object-by-object surveys, the construction of composite SEDs from galaxies 
spanning a range of redshifts allows for the cr eation of very high-quality and wel l-sampled 
SEDs that will be invaluable for SPS studies (|Assef et alj|2008l jKriek et al.ll201lh . 

Another theme of this review has been the growing realization that uncertainties in the 
SPS models are becoming a critical limiting factor to the interpretation of galaxy SEDs. The 
challenge here is not simply to enumerate the uncertainties but rather to identify areas where 
clear progress can be made. As a first step, all SPS models should include contributions 
from nebular emission and dust around AGB stars, as these processes are known to occur 
and the incorporation of such effects into the models is reasonably straightforward, even 
if the details are uncertain. Panchromatic models (i.e., FUV-FIR coverage) should also 
become standard both because IR data are no w widely available and because sophisticated 
dust emission models are well-developed (e.g.. iDraine fc Lill2007r i. The stellar atmospheric 
and synthetic spectral models will benefit from new asteroseismology measurements from 
the Kepler mission, interferometric observations (e.g., by CHARA), and new very high 



resolution UV-NIR spectra l atlases of nearby stars across the HR diagram (jBagnulo et al. 
20031 . iLebzelter et al.ll2012T ). Perhaps the most vexing issues lie with the stellar evolution 



uncertainties, as obvious calibrating data are lacking. The well-known problem is that, while 
globular clusters are the canonical testing ground for stellar evolution, metal-rich clusters 
are rare, and so the models tend to be poorly constrained precisely in the metallicity range 
most relevant for modeling galaxies. Moreover, the largest uncertainties are associated with 
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fast evolutionary phases and so stars in such phases will be rare in all but the most massive 
star clusters. Efforts to constrain uncertain stellar evolutionary phases from the SEDs of 
galaxies is promising because the right metallicity ranges are probed and there are sufficient 
numbers of stars to overcome Poisson noise, but the obvious com plexity of dealing with 
composite stellar populations makes t his approach challenging (e.g., MacArthur et al. 2010l . 



Kriek et al. | 201C j_ ZffjettietajJ I2012T ) . The Panchromatic Hubble Andromeda Treasury 
survey ( DertcantoneTlrtT 2012 ) is an HST program covering ~ 1/3 of M31's star-forming 
disk in six filters. It promises to provide new and powerful constraints on luminous and 
advanced stellar evolutionary phases at moderately high metallicities. 

More accurate models and higher-quality data will necessitate a more sophisticated ap- 
proach to comparing the two. Presently, model fitting is something of an art, owing to 
the fact that large regions of parameter space are often severely under-constrained, which 
implies that the choice of priors on the model parameters can have a large impact on the 
derived results. A few basic guidelines should be followed to ensure that results are ro- 
bust. For example, one should not simply fix a parameter to a particular value because 
it is under-constrained. Moreover, because the likelihood surface often contains multiple 
peaks and valleys and is frequently computed on a coar se grid, the best-fit parameters ought 



not be chosen based on the minimum of \ (see e.g.. lTavlor et al.ll201lf ). Rather, the full 
posterior distributions should be used to derive best-fit values and associated uncertain- 
ties. The choice of priors needs to be considered carefully, and in fact the model space 
should probably depend on the type of data being fit, the redshift of the object, and even 
its spectral type (quiescent vs. star-forming vs. peculiar). As the number of parameters 
increases, Markov Chain Monte Carlo techniques will see more widespread use owing to 
their efficient exploration of parameter space. With regards to the modeling of moderate 
resolution spectra, the analysis of spectral indices should eventually give way to full spectral 
fitting as the latter not only allows for the extraction of more information but also allows 
the modeler to visually inspect the fits in a way that is not possible when only EWs are 
extracted from the data. This is important for identifying model systematics and areas for 
future improvement. 

Finally, further work is needed to understand what is knowable, in principle, from the 
modeling of galaxy SEDs. Questions such as 'how many discrete SF episodes can be mea- 
sured from high quality optical spectra?', or 'how many moments of the metallicity distribu- 
tion function can be extracted from SEDs?', or 'can the detailed dust attenuation curve be 
measured with sufficiently high quality data on an object-by-object basis?' have yet to be 
thoroughly explored. Addressing these questions will be difficult because they depend sen- 
sitively on the quality of the data, the SED type, and the reliability of the models. Fitting 
routines such as STARLIGHT, STECMAP, and VESPA that attempt a non-parametric 
recovery of the SFH and metallicity distribution function offer probably the most reliable 
tools to explore these questions. Theoretical studies aimed at understanding 'what is know- 
able' will help guide the next generation of surveys aimed at studying the detailed physical 
properties of galaxies. 
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